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1.0 INTRODUCTION 


This manual is written for users of tha threo-di mensi onal version of tha TEACH computer 
program as developed at P&UIA. 

This program solves the governing transport equations in Reynolds average form for the 
flow of a three-dimensional, steady state, viscous, heat conducting, multiple species, 
single phase, Newtonian fluid with combustion. The governing partial differential 
equations are solved in physical variables in either a Cartesian or cylindrical coordinate 
system. The effects of rotation on the momentum and enthalpy calculations are modeled in 
cartesian coordinates. The flow of the fluid should be confined and subsonic with a maxi- 
mum Mach number no larger than 0,5. 

The program is flexible and is written in a modular form. By simply turning the calcu- 
lation of various transport equations in the input on or off, several types of flow can be 
solved. For example; either laminar or turbulent, reacting or nonreacting, single or mul- 
tiple specie, and isothermal or noni sothermal flows can be solved. The fluid may contain 
up to four distinct species, which may be inert or undergoing turbulent combustion. 

In order to demonstrate the range of fluid flow problems that can be solved by the above 
procedure, several sample cases are included. These cases range from simple laminar pipe 
flow to the flow field generated by a JT9D cooling louver. 

The program has been made quite versatile by adding a useful feature which allows blockage 
inside the flow field to be defined by the input. This feature enables the flow inside 
almost any geometry, within the constraints of the coordinate system, to be computed. 

An effort has been made to include ’warnings' etc., so that the solution obtained is the 
numerically correct representat i on of the flow field. However, Computational Fluid Me- 
chanics is, in spite of the advances made, a black art and the user has to rely on his ex* 
perience to ensure that an accurate solution has been obtained. The following criteria 
should be helpful in this respect. 

1) The input is properly defined. Close inspection of geometry plots and careful examina- 
tion of every variable field in the printout is necessary to make this determination. 

2) The solution has converged. All the residuals are smaller than the error determined to 
be acceptable for the case. Determination of an acceptable error is usually done from ex- 
perience, and generally varies from 0.5& to 5 X, 

3) The solution has stabilized. Flow properties in sensitive regions of the computat i onal 
domain have stopped changing from one iteration to another. This criterion, in conjunc- 
tion with the convegence criterion, is necessary to ensure that a numerically correct sol- 
ution has been obtained. However, in some cases, where the convergence criterion is 
difficult to satisfy, acute monitoring of the variable fields during the iteration process 
can be relied on for an acceptable solution. 

4) A grid independent solution is obtained. Grid i ndependence is attained when the use of 
a finer mesh does not change the results significantly. In a 3D code, the number of nodes 
available is often not adequate to go to a finer mesh. In such situations, some idea of 
mesh independence can be obtained by going to a coarser mesh. 

5) Mesh_Reynolds numbers and Peclet numbers are within limits. The wall Reynolds numbers 
and Peclet numbers are printed in the output. Wall Reynolds numbers should be between 30 
and 400. Peclet numbers should be between -5 and +5 in regions where any field variable 
has a high gradient. 

6) Cell aspect ratio is acceptable. The aspect ratio of computat i onal cells should be 
unity to ensure accuracy, rapid convergence and to avoid numerical instability. However, 
in most cases this criterion may not be satisfied and under relaxation factors may have to 
be reduced significantly to achieve convergence. 

Since proper definition of the computational grid is essential for smooth convergence and 
accurate solution of the flow field, the user is advised to make a preliminary run with a 
coarse grid. The flow field obtained should then be inspected carefully to ensure that 
the boundary conditions are properly specified and to ascertain regions of steep gradients 
which will require a denser grid. The problem should then be run again with a fine grid 
using a denser grid in the regions of steep gradients. An option is provided for the in- 
spection of the grid before running the program to facilitate the selection of optimum 
grid spacing. 
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In the initial running of TEACH, it is often difficult to obtain stability in the simu- 
lation of turbulent reacting flow. Guidelines to aid the use are given in Appendix F. 

The program can bo run interactively or submitted to batch from any CMS terminal. Running 
interactively allows the user to check the geometry and boundary conditions and make a 
short run to see if a case will start properly and examine intermediate calculations to 
determine if it is going to solve the problem intended. Because it is not always practi- 
cal for the user to run all cases to completion on CMS, a version of TEACH is available to 
execute these cases as a batch job. The decision to run batch or CMS is made at the be- 
ginning of a CMS terminal session. 

The user will require some basic experience before he is able to employ all the capabili- 
ties of the program. It is strongly advised that a few simple problems be solved before 
attempting to soIvq a complicated problem. In this regard reproducing at least a few of 
the samplQ cases solved in this manual might prove useful. 

In the next section the theoretical foundations of the computer program, which include the 
governing transport equations, the discretized equations, the solution procedure, and mod- 
els of combustion and turbulence, are described briefly. See reference 1 for details. 

In Section 3 the mechanical structure and operation of the program are described. This 
includes a functional diagram, a flow diagram, line by line description of the input and a 
sample of the output options available. A typical CMS session is also included. 

The final section contains the sample cases. Each case includes a copy of the input, grid 
plot, sample output, and streak line plots in the planes of interest. 

2.0 ENGINEERING DESCRIPTION 

This section presents the theoretical foundations of TEACH to familiarize potential users 
with the TEACH model, and includes the coordinate system, the basic differential equations 
being solved, and the finite difference versions of those equations used in the model. 

2 . 1 Grid and Coordinate System 

TEACH uses a 3-dimensional orthogonal coordinate system (Cartesian or cyl i ndr i cal ) . The 
grid formed by the intersection of coordinate lines is used in the modeling of fluid flow 
systems. The intersection of the grid lines form the grid nodes at which all flow proper- 
ties except the velocities u, v, and w are calculated. The axial velocity, u, is calcu- 
lated midway between grid nodes in the axial direction, the radial velocity, v, is 
calculated midway between grid nodes in the radial direction, and the tangential velocity, 
w, is calculated midway between grid nodes in the tangential direction (figures 1 and 2). 
Directions in the grid are identified as north, south, east, west, front and back. 

A breakup results in defining an array of rectangular cells, each bounded to its west and 
east by X-gridlines, to its north and south by Y-gridlines and to its front and back by 
Z-gridlines. The u, y, and w velocities lie in such a manner that normal velocities lie 
on scalar cell boundaries (Figure 3). 
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2.2 Governing Differential Equations 

Solution of a flow problem in TEACH consists of solving a set of equations including# the 
continuity equation (conservation of mass)# the three equations of motion (conservation of 
momentum)# the energy equation# the thermodynami c equation of state (relating pressure# 
density# and temperature) # and a suitable set of scalar equations such as the species 
transport equations. These equations are derived for nonreacting# compressible# multiple 
species# laminar flow. For an explanation of the turbulence model# the combustion model# 
and their inclusion in the equations developed here# see reference (1). 

Laminar flow problems can be solved with the Navi er-Stokes equations and a suitable set of 
scalar transport equations. The transport equation in its general form states that the 
net amount of a scalar such as species or energy carried into a control volume by the mass 
entering it (the convection term) plus the amount entering through diffusion or conduction 
(the diffusion term) plus the amount generated within the control volume# such as through 
a chemical reaction (the source term)# must equal zero. Turbulent flow equations are sim- 
ilar to the laminar flow equations except that viscosity is replaced by an effective 
viscosity which includes the "turbulent" viscosity caused by the turbulent eddies# and an 
effective pressure term due to the turbulence kinetic energy. 

2.2.1 Ihfl_ Cent inultv Equation 

The continuity equation is a special case of the transport equation in which the mass it- 
self is being balanced# there is no diffusion term and no generation term. The continuity 
equation expresses the fact that for any control volume under steady state conditions the 
mass flowing into the volume must equal the mass flowing out. The flow field is specified 
by the velocity vector 


V * 1u + jv + kw 


where i# j and k are the unit vectors# and u# v# and w are the components of velocity in 
the X-# Y-# and Z-directions respect i vely . 

The mass flowing out of the two faces of the control volume p o erpendi cular to the X-axis 
having area dydz are given by 


(2) pudydz and (pu + dx) dydz 


Similiar expressions can be written for the front and the back faces. Subtracting the 
mass flowing into one face from the mass flowing out of the other gives the net mass flow- 
ing into tho element in the X-direction 


(3) dxdydz 


A similar expression can be written for the mass flowing into the element in the Y-# and 
Z-di rect i on . 


|y“ dydxdz 


(5) 


dzdxdy 
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Addition of equation 3, 4 and 5 givos the equation of continuity which states that the 
mass flowing into the volume must equal the mass flowing out for a steady state system 


C6 ) 


3 P° + 3pV . 3pW _ n 

ax 3y Jz~ ~ U 


2,2,2 Momentum Transport Equations 


The equations of motion are derived from Newton f s Second Law* which states that the summa- 
tion of forces acting on a body is equal to the time rate of change of its momentum (i.e. 
mass times velocity) 


(7) 



d(mv) 

~ST~ 


considering the right hand side of the equation first# the mass of fluid entering the west 
face is P udydz. The velocity of this fluid is u and therefore the momentum per unit time 
entering the west face is 

o 

(5) (pudydz)u = pu dydz 



FIGURE 4 

U-MOMENTUM TRANSPORT INTO A FLUID CELL 
(For clarity/ not all flows are shown) 

Figure 4 shows the expressions for the U-momentum crossing the faces of the Control 
Volume. 
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Fluid flows into the control volume through the south faco at the rate vdxdz with a ve- 
locity component in the X-direction of u* so the contribution to the X-momentum per unit 
volume is 


(9) (pvdxdz)u 3 puvdxdz 


The momentum per unit time leaving the east face is 


(10) [pu 2 + (pu 2 ) dx] dydz 


Similarly* the momentum per unit time leaving through the north face is 

(id [puv + -i (puv) dy] dxdz 


Similar expressions can be written for the front and the back faces. Subtracting the mo- 
mentum per unit time entering from that leaving* gives the net increase in X-momentum per 
unit time of the fluid in the control volume 


(12) [ ■§- ( P u 2 ) dx] dydz + [ \rj (puv) dy] dxdz + [ \= (puw) dz].dxdy 


Any increase jn momentum of the fluid in the X-direction is produced by forces acting in 
the X-direction on the control volume. In general there are two types of forces to be 
considered; Body forces such as gravity and surface forces such as pressure and viscous 
shear . 

When a control volume having dimensions dx* dy* and dz is isolated from a fluid* the forc- 
es acting on it while in the fluid can bo represented as shown in figure 5. 


t°yy + dy] dxdz 



(For clarity* not all forces are shown) 
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Gravity is a body force acting on every particle in the control volume. Its contribution 
to the external forces acting on the control volume is 


(13) 9 x pdxdydz 


where g is a vector denoting the gravi tat i onal force per unit volume. If the coordinate 
system is rotating, centrifugal and Coriolis forces must be added as body forces. 


(13.1) p [ 2 uxv + ux(uxr) ] 


Where u i s the vector representing the angular velocity of rotations, » is the vector re- 
presenting the velocity of the fluid and fis the vector representing the location of the 
center of the fluid cell. 

U = T<*) x + jUy + 

7=1 u + 3v + 
r = Tx + 3y + 


See reference (4) for the development of (13.1). The other components shown in the figure 
are pressure and shear forces which act only on the surface of the control volume. The 
x-momentum enters and leaves the control volume as a result of normal stress <7 X x and 
shear stressdyx andd . Those stresses can be expressed in terms of strain rates in the 
following manner. zx 


a =-P + 2 M iH +x (£ii + ^. + «) 
xx M ax ' 3x ay a z 


au . av A aw 


. / av , au i 

U5 > °yx " » ( 17 + ly > 


and 


(16) 


'ZX 


m ( iSt + i! i ) 

M ax az 


where M is the viscosity of the fluid and A. is the second coefficient of viscosity which 
has an empirically determined value of -(2/3) 


The x-momentum entering and leaving the control volume by the action of these stresses can 
be treated in a manner similar to the x-momentum entering and leaving the control volume 
by convection. Hence summing all the forces in the X-direction and dropping terms that 
cancel yields 


(17) 


( v 


3o 


XX 


-ax 


a a, 


£X 


3a 


ZX 


ay 


az 


) dxdydz 


Putting equations 14, 15, and 16 in equation 17 yields 


| 9 X P " t£ + h- C 2 m 


3X 


ax 


au 

ax 


o ax ay az /J ay ay ax n 


av 


aw 


az 


dxdydz 


(18) 
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Equating the rate of change of momentum to the forces acting on the control volume C ex- 
pressions 12 and 18 respectively) and dividing throughout by dx dy dz we have the follow- 
ing x-di recti on equation of motion 


(19) 


a(puu) + a(puv) + 
ax ay 

+ i_ ( M ( iE + jv 

ay ay ax 


a( puw) _ aP + a 

az ” P9 x ” ax. ax 




jW + _3U M 
3X dZ 


ro 3U 2 f 

[2 “ 77 - 7 1 


du_ + 11 + 2* ) i 
3X ay a z J 


2*2,3 Thermodynamic Equation of State 


Density at every point in the control volume is calculated using the equation of state 



where hi | is the mass fraction and Myi is the molecular weight of species l . 


2.2.4 Scalar Transport Equations 

An example of a scalar transport equation is the energy equation which has a convection 
term, a diffusion term, and heat generation terms. The convection term represents the net 
heat carried into tha control volume by the mass flowing across its boundaries at a rate 
proportional to the velocity normal to the boundary. 


<21> -fx ^P UH ) + ly (pVH) + TZ (PWH) 


The diffusion term represents tha net heat entering the control volume through conduction 
at a rate proportional to the thermal conductivity of the fluid 


C22> If <K £ ) + h (K -£> ♦ «•!> 


ay 


ay 


32 ' az' 


and heat carried in by mass diffusion 


(23) 






The generation term S(H) consists of radiation and shear work, and when the corrdinate 
system is rotating, S(R) consists of work done on the fluid by centrifugal and Coriolis 
forces. The terms due to centrifugal and Coriolis forcas are: 


(23 5) S(R) = E w ( cj . r) - r ( cj • u ) ] 


See reference (4) for the development of S(R). See reference (1) for the development of 
the shear work terms and reference (2) for development of the radiation terms. The final 
energy equations becomes: 


if (ouH1 + 1/ (,>vH) + 71 ■ h «77> + 7y < > * h <4> 

+ |y t |_ ( 2ri hj, |2i ) + S(H) t SIR) 



( 24 ) 
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The specie transport equation must be solved for each species present in the control vol 
ume. The convection term represents the mass fraction flowing into a control volume 


(25) Tx ( pura l) + !y + Jz {pwm i* 


The diffusion term represents molecular diffusion of the species into the control 


volume 


C 26 ) 




where 


Tg is the mass 


diffusion coefficient of species 1. 


The source term 


(27) (T, mj) 


represents the time averaged reaction 
ature and mass fraction of the species. 


rate of species mg which is a function of temper- 
The final species equation becomes: 


( 28 ) („mjt> ♦ i- fpvmjt) ♦ Um t ) * I* iSl ♦ Jy r< BL ♦ M * S £ (T. , t ) 

2.2.5 The Modeled Equations 

To summarize the equations being solved and the flexibility afforded by the program, as- 
sume that a very basic problem is to be solved. Such as a constant density, constant tem- 
perature, single species liquid flow field. The unknowns in this example are u, v, w, and 
p. The three equations needed are: 

the continuity equation 


(29) 


3X (pu) + ly (pv > + It (pw) 


az 
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the X-momentum equation 

a(puu) + a(puv) + a( P uw) 


ax ay 32 

(30) + i_ ( u ( l!i + iv )) + i_ („/ 3W . 3U n 

3y ay ax lz lM( lx TI 


«n 3P . a /0 au , 

P9 x ax + ax ' 2p ax ) 


+ p[(v - Z u y + Xu 7 ) u - (w - Xu + y u ) u + Vo - W u ] 

* y * y z y 


X 2 

the Y-momentum equation is 


jlpuv) * 3( pVV) * 3( c vw) »P . 3 , , 3U . 3V ,, 

SX 9y 3z " 9 y 17 « ( “ 1 '3y + u )) 


C 31 ) 


+ 4^ 44)) 


ay 


ay 


a z’ 


ay 


az 


+ p[(w - xuy + yo x ) o x - (u - yu z + zUy) o z + wu x - uu z ] 


the z-momentum equation is 


a(puw) + a(pvw) + a(pww) 
ax ay az 


„ aP + 3 („t aw + 3,1 n + i- lul — + — )) 
pg z " Iz + lx (ll( lx + lz n + ly Ul az ay 11 


3 V . 3W 


(32) 


+ i_ (2 M 1* ) 

32 1 * az 


♦ P [(u - yu z + zu y ) Uy - {V - zu x + XU Z ) u x + Uu y - Vc x ] 

If the problem is changed to a non isothermal gaseous flow extra unknowns, P and T are to 
be calculated. The equation of state and the energy equation must also be solved. 


(33) p . p RT 


Energy equation 


(34) 3X 


(puH) + i_ ( pvH) + -L. ( pW H) 


32 


h < K £> *Tf «£> +h u Ti > + S 


If the flow field involves different specie undergoing mixing, the energy equation has 
more terms added to account for the heat carried by the diffusion of the individual 
speci e. 


Energy equation 


jj- <PU ») + <p vH > (pwH) 


(35) 


ay 


32 


3 

3X 


(Kil) 


+ i_ (k il ) 
ay ay ' 


+ — (K 
az 


_aT 

az 


+ |f ) + *= (zr^lf ) ♦ £ |f ) * s 


where a species transport equation must be solved for each species, 1. 
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(36) i_ (pumj) + (pvmjj) «■ (p wm i) - h (r i 77^ 1 + Ty (r i ) * 77 ( Q if } 

For turbulent flow# two. more equations must be solved; one for K , the turbulence energy , and 
one for € , the energy dissipation. The turbulence model replaces the laminar viscosity 
MI with the eddy vi scosi ty Meff i n the momentum equations. It also adds extra terms to the mo- 
mentum equations, the species transport equations, and the energy equation. For reacting 
flows, a reaction rate term which determines the rate at which the fuels and oxidants react 
to form products is added to the specie transport equations, equation 36, as a source term. 

All of the transport equations can be algebrai cally manipulated into the form: 

<37, (pu „, ♦ l. („<) * i. M) ■ r j£ * ^ ri£ ♦ r j| * s, 

where the left hand side contains the convection terms and the right hand side comprises 
the diffusion terms and the source terms (5). Any term that is not a diffusion term or a 
convection term is added to the source terms. Table 1 contains a tabulation of the terms 

in r and $ . 


i r s. 


Conti nuity 
x- momentum 

y-momentum 

z -momentum 

Energy 


1 0 0 


u u 


v u 


H K 


3 3U 2 3U ^ t 3W + 2_ 3V 4 . 3 3W 3P + 


W U TTf TT- 


3 

3U 

2 au 

+ 21 

+ 21 

3X U 

3 X 

7 u 37 

ay 

az 

+ P 

C(v - 


xu z ) 

°z " 

3 

3 V 

2 au 

. av 

. aw 

Ty U 

ay ' 

7 u 77 

ay 

+ -- 

3Z 

+ P 

• — i 
* 

Xu y + 

y “x* 

u x - 1 

3 

3W 

2 au 

+ av 

+ 21 

3Z M 

3 Z 

7 M ax 

ay 

3 Z 

+ P 

C(u - 

C 

1 

+ 

Zu ) 

(J - 1 



£ 

y 

y 


ay 


3X 3Z 


3X 


(w - xu y + yu x ) u y + v« z 


3 

3X 


3U 

ay 


3 

3Z 


aw 

ay 


3X p 9x 
Wu y ] 

a P . 

3v P ^y 


yu z + Zu y ) u z + Wu x - Uu z ] 


a 

3 X 


3U 

32 


3 

ay 


3 v 

3 Z 


3P 
3 Z 


(y - zu x + xc> 2 ) u> x ♦ Uuy - y<*> x ] 


3 /h am\. 3 (h 3m \ + 3 (h 3m \+ s 
3 X y Tx) + 3y ( h Ty J Tz TTj S T 


- P |u [u y (yu x - Xu y ) ♦ u z (zu x - xu z )] ♦ V [u x (Xu y - yu x ) 


u z (zu y " yu z n + w tu x (XM z 


z« x ) + u y ^ u 2 


Zu 


,)]} 


Specie transport r* 


Table 1 
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2.3 Discretizing the Equations 


This section contains a presentation of the technique used by TEACH to solve the general 
transport equations, given in section 2.2.5. The goal is to describe the approxi mat i ons 
that are made in derving the finite difference form of the equations from the differential 
equations. This knowledge will allow the user to recognize the errors that could arise 
and enable him to avoid them. There are several techniques that can be used to discretize 
the differential equations all of which produce essentially the same difference equations 
(see reference 3). The technique used here is the integral method. In this method the 
differential equation is integrated over the volume of a cell. The cell geometry given 
below will be used throughout this section. 


FIGURE 6 

SUBSCRIPT NOTATION FOR A FLUID CELL 
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Integrating the convection terms over the cell volume: 


(38) 


y i elds 


(39) 


separat i ng 


(40) 


f z b r y n /* x e 

\ r *y X 3 X dxdydz 

• » w 


/ f { (pugf) _ ( pu ^ 1 dydz 
z f y s e * w I 


f h /*" f Zb / y n 

-4 J. (pu ^x d * dz -/ J (pu <H dydz 


Z f ’ y 


6 ' z f y s 


At this point Me make our first approximation; P, u, and $ are all functions of x, y, and 
Zt but using the mean value theorem, which states 


(41) 

whero 

(42) 

y i elds 

(43) 
but 

(44) 


/ 


A2 


f(z)dz * f ( z) AZ 


AZ 


z--|i<z<z+|^ 


( pUflf ) 


x e, y p z p ’ y s } (zb " zf) " (pU<J) x w p p z p (y n “ y s ) * (zb “ zf ) 

(y n - y e h* ( z k - a and - yj* (z. - zJ - a 


n J s'x a b ‘f # x a e 
e e 


n J s'x w b f x °w 
w w 


where a e and a w are the surface areas of the "east and west cell walls respect? vely, 
V i eldi ng 


(pU<0 x y z a e 
x e* y p’ z p e 


{f>yi) x y z a w 
w’ y p ,z p w 


(45) 

Similarly the other convection term works out to be 

(46) 


(pvd) _ a n 
x p ,y n’ z p n 


( pv«0 a 

x p’ y s’ z p s 


(47) (pwd) a. 

V y p* z b b 


( pW <5) a f 

p ,y p’ z f r 
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Integrating the diffusion terms over the cell volume 


(48) 


yi elds 


(49) 


P /" p *_ r * 


3X r 3x dXd > dz 


w 


C C PP«¥. 


) dydz 


Using the mean value theorem again gives 
<50> 

by using the midpoint rule (see reference 3) we can approximate d$jdx and get 


(51) 


’ V 

Ay 


"n *n 


Up - W 

Ay 


r s a s 


Si mi larly 


(52) 



AX 


r e a e 


Op - , 

AX S# \ 


and 


( 53 ) 


<"b - V „ Op - *F> . . 

n b a b 4Z f *f 


1 1/16/84 

03/04/86 
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The source term is linearized as: 

Srf ■ % + B 


where A and B are assumed contstant. Integrating the linearized source term over the con 
trol volume as* 




(Ad p +B) dxdydz = Ad p . vol + B 


vol . 


Def i ni ng 


A Vol =Sp and B 


Vol=Su we have 




dxdydz = s p d p + s u 


The transport equation can now be written out 


C54) ” C lAr + C N**n " C S^s + C B^B " C F*V = °E^E " ^ " D W^P " 

+ D N (r5 N - 6p) - D s (d p - d $ ) + Dg (^g - dp) - Dp (dp - dp) + Sjj + Spdp 

where Cg* C|sj» etc, are the convection coefficients 


(55) c E - ( P u) e a e , C N * ( P u) n a n , etc. 


De f / etc, are the diffusion coefficients 


r r 

(56) D r * — a , D m * — a , etc. 
E ax e ’ N Ay n * 


and SU comprises all the source terms which do not contain 4>p and Sp<?>p comprises 
all the source terms containing <J>p . 
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Noto that in equation 54, 0 w , 0 e , 0 f » 0 , 0 s , and ^ n are unknown * 



.*N 







*w. 

• 

4> p 

4>« * 4>z 



* s 




*s 


.*» 






FIGURE 7 

LOCATION OF 0s USED IN TRANSPORT EQUATIONS 

0 is assumed to vary according to tho following relation between cell centers, so the un- 
known 

0 5 can be calculated from the known 0 s according to 


(57) * g E + (l - g E ) d p 


whero g is a weighting factor defined as: 
g = 1/2 - 2 1 Pe $ 2. 

g = 0 Pe >2 

g = 1 Pe <-2, where Pe = CE/DE 

The transport equation can be rewritten now 


C E ( Vp ♦ (1 - W - Wu + (1 - W * Wp * 11 - f n ,)<l N 
- c s (f s <( s ♦ (1 - f s )«t p ) ♦ C B ( f b rf p + (1 - f b)<B ) - Cp( f f «* B - (1 - ff>V ■ Vp * s u 


where 
(59) f 

e 




with algebraic manipulation and subtracting 0 times the continuity equation 


( C E ‘ C W + 


C N - 


C S + C B - 


C F ) ~ 


(60) 
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the transport equation reduces to 

~ ffp ) “ ^p) Ap(^p “ ^p) ” ^p) Ap(dp ~ tfp) 

+ Sp*(p + Sy ■ 0 

where 

“*> A„ « C„ (f n - 1) , A $ ■ C $ f s , etc. 

which can bo further simplified to the finite difference equation for 0: 

C63) Vp “ Vn + Vs + Ve + + Vb + + S U 

where 

(64) Ap = A M + A $ + A„ + A E + Ag + Ap - S p 


To soIvq the above equation, the weighting factors fe, fw, fn, fs, fb, and ff have to be 
calculated based on the values calculated for g , g n , g s , gh, gf, and g w for each cell. 
The method. used to calculate weighting factors is given in reference 1. 


2.4.1 


If the cells are numbered from 1 to HI in the X-direction, from 1 to NJ i n the 
Y-direction, and from 1 to NK in the Z-direction, the equations will be solved in the 
field from 2 to HI— 1 » from 2 to NJ-1, and from 2 to NK-1. The values lying on the border 
of the domain are set by the input boundary conditions. There are three types of boundary 
conditions which can be applied. One type is where the value of <p is specified, such as 
the velocity in an inlet. Another is where the gradient of <p i s zero, such as the zero 
temperature gradient at an adiabatic wall. The last is where <t> is unspecified such as at 
an outlet. On an axis of symmetry, v-velocity and the gradients of all $5 except 
v-velocity are set to zero. The values of field variables in the boundary cells adjacent 
to a periodic boundary are adjusted every iteration by making them equal to the corre- 
sponding cell next to the periodic boundary on the opposite end of the control volume. 
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The handling of wall effects on heat transfer, turbulence energy, and energy dissipation 
rate is beyond the scope of thia manual. Refer to reference 1 for details. 

2.4.2 Description of the Solution Algorithm 


The finite difference equation to be solved is the transport equation 

(65) Apgf p = A h 6 h + A s d s + A w jJ w + A E d E + Agrfg + Aprfp + Sg 

where can be u, v, w, P, k, e. Mg, or H. The equation gives the value of (|) at a point P 
as a function of its neighbors to the north, south, east, west, front, and back and must 
be solved for each cell in the flow field. The method used to solve all the equations is 
to solve for values of (J in a given plane one row or one column at a time. If the fields 
are solved row by row then the values in the neighboring columns and planes are assumed as 
constant. Once one plane is solved, then the neighboring plane is solved until the whole 
block matrix is updated. 


Vp ■ Vm + Vs + c ' 

or 

(67) Vj = Vj+1 + Vj-1 + c ' 

where 

(68) C - A^ + A e «( e + A B d B + Apdp + Sy 

In this case all the 4>s in one column are calculated at once using a Tri-Diagonal Matrix 
Algorithm (TDMA), and solution of the whole field proceeds by solving each column sequen- 
tially in this manner. 0 

2.4.3 Under relaxat i on 


Because the finite difference equations being solved are all non-linear, allowing the pro- 
gram to calculate and assign new values of <f> in one set of finite difference equations can 
have a drastic effect on other variables being solved, leading to convergence problems. 
To prevent this, the values of <J> calculated are under-relaxed before they are used in oth- 
er equations. 


(69) 



calculated 


(l-F)rfp 

old 


inhere F is an under-relaxat i on factor usually between .3 and .7 depending on the variable 
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2.4.4 Pressure Correction 

At this point it should be noted that none of the transport equations can be used to cal- 
culate pressure, yet the momentum equations use it. Another equation is needed to solve 
for pressure. The extra equation added is the pressure correction equation and the algo- 
rithm setting it up and using it is called SIMPLE (Semi Implicit Method for Pressure 
Linked Equations). Briefly, this method does not calculate a pressure field but corrects 
a pressure field from an initial guess or from the previous iteration. 


At the start of an analysis, the pressure field is assumed to be constant and the velocity 
fields are solved using the momentum transport equations. Because the pressures used in 
the momentum equations are guesses, there will be errors in the velocities calculated and 
continuity is not conserved. Using the new velocity fields the continuity equation is ap- 
plied to each cell and the error in the net mass flow into the cell is saved as a "source" 
term. 


(70) Mp * (pu*) e « e 


( pu *) w a w + ( pv*) n a„ - <pv*) s « s + «b 


( pW* ) f a^r 


If us subtract this from tho continuity equation 

(71) (pu) e a e - (PU) W + (pV) n a R - ( pv) g o $ + (p«) b - (p«) f a f * 0 


we get 


(72) -Mp ■ (pU') e a fi - (pu') w + (pV') n a n - (pV') $ o $ + ( pW* - (pw')f 
where 


(73) 


u f is the velocity error, u is the exact velocity, and 
the solution of the momentum equation. 

Similarly v f and w 1 can also be found. 

If we assume that the velocity errors between two cells 
ference in pressure errors between those cells, we can 
constant of proport i onal i ty for the velocities in the 
pressure correction equation: 

(74) ApPp - A n P n + A $ P S + + A e P e + AgPg + ApPp + Sy 

where 

(75) A P “ A N + A S + \ + A E + A B + V 


A n f A s , A w, A e , Af, and A^ are the constants of proport i onal i ty which come from the momen- 
tum equations, and 


U* is the velocity obtained from 


is proportional only to the dif- 
substitute those pressures and a 
continuity equation and get the 


(76) 



The pressure correction equations are solved and the pressure corrections are added to the 
pressure field. The velocity fields are also corrected using the assumed proport i onal i ty 
between pressure difference and velocity. For more information about SIMPLE see reference 
1. 
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2.4.5 Convergence Criteria 


Convergence is determined at the end of each iteration on the basis of the 'residual 
source' criterion. The residual source R^ is defined as 


(77) = A p**p ~ a n^n “ A s^s " a e^e " " a b^b " a f^f " s u 


For convergence# R 


(78) 


ZM < p 

■R < e 

Vef 


summed over the entire flow field: 


where R<f>ref is a reference flux based on the total fixed flux of $ into the domain and the 
error is of the order of .001. 


When mixing calculations are involved an additional requirement that the sum of the mass 
fractions at each node be less than (li 6) must be satisfied. When variable density flows 
are being calculated there is an additional requirement that the change in density at each 
node in one iteration must be less than e. 

2.4.6 Solution Procedure 


1. Guess fields for all variables. 

2. Solve for u* # v*# and w* using prevailing pressure field. 

3. Solve pressure correction equation and correct pressure and velocity fields. 

4. Solve equations for all other variables. 

5. Update fluid properties such as density and viscosity. 

6. Test for convergence. If not converged# use prevailing fields as guesses and repeat 
from (2). 
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3.0 PRPQRAfl PE$CfU PT IPN 


3 . 1 Functional Diagram 


This section contains a functional description of each major subroutine in the program. 
Following each subroutine name (given in capitals) is a brief description of the subrou- 
tine's purpose in the program. Branches emanating from the subroutine name denote rou- 
tines which are called by that subroutine. 


MAIN- 


-CONVER- 

— GEONTY Expand block geometry input to cell geometry. 


(-INIT- 


KRESTRT- 


Controls program execution 


T Read in ar 

NABOP 

5ETDFT 

XINPUT 


(—INPUT i Read in and print out all input cards. Set defaults 

Check boundary conditions of fluid blocks 
Set all defaults for the case 
Explain the input in English 

Convert units from English or SI units to engineering units* sets program constants. 


(-CELLST- 


Initialize all boundary conditions. 

Initialize the flow fields if not a restart case. 


— Set up boundary condition tables. 

SETUP Set up table for each boundary condition type 

■ Initialize enthalpy field from temperature field. 

(— ENTHMX - Calculate enthalpy of mixtures when en 


(— HFROMT 


table is used 
Store all field values in restart file. 


enthalphy 


— r — store ail t 
I— GETPAM 


Get restart file from PANVALET 


•GRIDfOUTLIN— Display expanded cell breakup, 
h PROPS 


|— PRINTR- 


— Calculate fluid properties 

(viscosity, density, specie diffusivity, specie emmisivity). 

— Print all variable fields. 

I— CONVER Convert engineering units to English or SI units for output. 

(— PRINT1 »PRINT2 Print initial variable fields 

PRINTH Print initial heat balance information. 


— CALCU Calculate finite difference equation coefficients 

— CALCV Calculate finite difference equation coefficients 

— CALCW Calculate finite difference equation coefficients 

— CALCP Calculate finite difference equation coefficients 

— DIVRGE Calculate divergence of velocity. 

— CALCTE Calculate finite difference equation coefficients 


for U-velocity field.# 

for V-velocity field.# 
for 14-velocity field.# 

for pressure correction field.# 
for turbulence energy field.# 


FIGURE 8 

TEACH Functional Diagram 


* Refer to Figure 8a 
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-CALCEO 

-CALCF- 


Calculate finite difference equation coefficients for energy dissipation field.* 
— Calculate finite difference equation coefficients for species mass fraction fields. 

-FSORCE Calculate extra source terms in the species mass fractions equations. 

-MODF Modify coefficients for species mass fraction boundary conditions. 

-TOMAIJ Solve the field using tri diagonal matrix algorithm. 


-CALCH- 


— Calculate finite difference equation coefficients for enthalpy field. 

-HSORCE Calculate extra source terms in the enthalpy equations. 

-MOOT Modify coefficients for enthalpy boundary conditions. 

-TOMAIJ Solve the field using tridiagonal matrix algorithm. 


-TFROMH- 


Calculate temperature field from enthalpy field. 

E NTH MX Calculate enthalpy of mixtures when enthalpy 

table is used. 


-CALCFH 

-CALCFL 


-CALCRD 


-PROPS- 


Calculate finite difference equation coefficients for mixture fraction field.* 

Calculate finite difference equation coefficients for mixture fraction fluctuation field. 
MOOFL— — — Modify coefficients for mixture fraction fluctuation boundary conditions. 
TOMAIJ— Solve the field using tridiagonal matrix algorithm. 

Calculate finite difference equation coefficients for radiation field. 

-M00RD Modify coefficients for radiation boundary conditions. 

TOMAI J * Solve the field using tridiagonal matrix algorithm. 


I -MCI 
-TOi 

l-MCN 
I -TOI 


-SEEIT- 


SEEIT2- 

-PRINTR- 


Calculate fluid properties 

(viscosity, density, specie diffusivity, specie 
Display convergence plot on Tektronix, 

Display convergence plot on A/N terminal. 


isiviiy ). 


-PR INTI , PRINT 2 Print variable fields 

-PRINTH— — Print heat balance information. 

-RE3TRT Save field variables in restart file. 

-SSOUTP User subroutine for special outputs. 

-PRINTN Print Peclet numbers and wall Reynold's numbers and wall shear stresses 

-FUR ATE Calculate reaction rates. 

-PLANES— | Select plane for streak line plot. 

STREAK plot streak lines. 


3—1 

1-STR 


FIGURE 8 Cconti nued) 
TEACH Functional Diagram 


x Refer to Figure 8a 
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Initialize any constants needed by this routine 

— MODX Adjust the coefficients along the boundaries 

— TDHAIJ ■ - Solve field variable using a tri-diagonal matrix solver 


1 1/16/84 

03 / 04/86 


FIGURE 8a (continued) 
TEACH Functional Diagram 
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3.2 Flow 


This section contains a flow diagram which shows the logical ordering of the functions de~ 
scribed in the previous section. A basic understanding of program flow will be helpful to 
users of TEACH in the interactive mode. 



FIGURE 9 

TEACH Flow Diagram 


















TEACH 3D 


PAGE 29 


ORIGINAL PAGE IS 
OF POOR QUALITY 


ISSUED: 11/16/8A 
REVISED: 03/04/86 


C BtSW ) 


Calculate Field Variables 


Calculate 
Prtifur* 
Coeff talents 


Adjust 

| eotffletanti for 
[ boundary conditio 


Solve 

Prtiim 

field 


is 

T 

Calculate 
0- vs loci ty 
Coefficients 

J 



Solve 

U- vs loci ty 
fisld 

D 


n 


Im 

T 

Celeulste 
V- vs loci ty 
Coefficients 

Li 

Adjust 



1 

p 


h 

boundary conditions 



T 

Celeulste 
H-vsloci ty 
Coefficients 





a 

fl 





ALTC 


Calculate 
Turbulence 
Coof f icients 

Calculate 
Coorgy dissipation 
Coeff talents 


Ad jut 

coefficients for 
boundary conditions 

I Adjust 

coefficients for 
| boundary conditio 


Sol vs 
Turbo Is 
fisld 


Sol vs 
Cnsrgy dissipation 
fisld 


d 

b 


Adjust U.V,snd U 
vs loci ti m 


F 
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3,3 

This section contains a detailed description of the card image input. Although the de- 
scription is complete it may be difficult to fully understand without an example to 
follow. Therefore referring to the sample geometry input in Appendix B may prove helpful. 

One of the most difficult decisions a user has to make while setting up the input to TEACH 
concerns the grid breakup. The user must trade numerical accuracy for CPU time, and ob- 
tain an acceptable solution in the minimum amount of CPU time. While the ultimate deci- 
sion can only be made after a process of trial and error, a few guidelines which apply to 
all cases are given here. 

CPU Time 

1) Solution time increases proport i onally to the square of the number of grid nodes. 

2) Given two similar cases with the same number of grid nodes but different grid spacings, 
the one with the most evenly spaced grid lines will converge sooner. 

Numerical Accuracy 

1) The most evenly spaced grid will give the most accurate solution. 

2) The finer grid spacings should be located near the walls, near the shear layers between 
coaxial jets, and near the interface between a jet and the surrounding fluid. 

3) Finer grid spacings should be located in regions where fluid variables have a high gra- 
dient, such as near an inlet and at a flame front. 

4) Mall Reynolds number (printed at the end of the output) should be between 30 and 400. 

5) Cell Reynolds number (printed at the end of the output) in any direction having high 
diffusion should be less than 5. 

For variables requiring units, either the English or SI system may be used. A summary is 
shown below. 



Engl i sh Uni ts 

SI Units 

Temperature 

°F 

°C 

Length 

ft 

m 

Veloci ty 

ft/sec 

m/sec 

Pressure 

lb/ i n 2 

N/m 2 

Turbulence Energy 

f t 2 /sec 2 

m 2 /sec 2 

Turbulence Energy 
Disipation Rate 

ft 2 /sec 3 

m 2 /sec 3 

Vi scosi ty 

lbm/(ft-5ec) 

kg/Cm-sec) 

Densi ty 

lbm/ft 3 

kg/m 3 

Enthalpy (Enthalpy File) 

BTU/lb-mole 

Cal/kg-mole 

Enthalpy (Output file) 

ft 2 /sec 2 

m 2 /sec 2 


The input to TEACH consists of one or two card image files. The main file comprises ever- 
ything required to solve a case if temperatures are not being calculated. If temperatures 
are being calculated and if specific heat is not constant, an enthalpy file must also be 
supplied. The enthalpy file has to be supplied if combustion is being calculated. The 
following is a description of the input for these two files. The enthalpy file is kept 
separate because one enthalpy file can be applicable to many cases. 
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3.3.1 Card Image Input 


TEACH uses a 3-dimensional orthogonal coordinate system (cartesian or cyl i ndr i cal ) . The 
grid formed by the intersection of coordinate lines is used in the modeling of fluid flow 
systems. Directions in the grid are identified as north# south# east# west# front and 
back# with the origin located at the southwest front corner unless an offset is input. 
With axi symmetry, an annulus is modeled by adding a radial offset to the Y coordinate so 
that the origin lies on the axis of symmetry# and the south wall of the breakup lies on 
the ID of the annulus. All blocks within the fluid volume must have the conditions on all 
boundaries (n# s# e# w# f# and b) described. 

Boundary conditions are input and stored as a table (see section 3.4). Each entry in the 
table is assigned a unique identifying integer label and a three letter mnemonic indicat- 


ing the 

type of boundary. A 

boundary can be one 

of six types: 


Boundary Conditions 

Mnemoni c 

Humeri c 

1 . 

Axis of symmetry 

SYM 

1 

2. 

Unspecified wall 

UWL 

2 

3. 

Specified wall 

SUL 

3 

4. 

Unspecified opening 

UOP (Nota 1) 

4 

5. 

Specified opening 

SOP 

5 

6. 

Periodic Boundary 

PER 

6 

Note 1: 

An unspecified opening can not be used 

in conjunction with an adjoining blockage 


Only one unspecified opening can be used in a setup. 

A specified boundary is one for which the user supplies all boundary values such as veloc- 
ities# temperatures# etc. An unspecified boundary is one for which the boundary values 
are calculated by the program. The specified boundary can be a wall or an opening. High 
shear forces are present near a wall# but not near an opening. To account for these shear 
forces# a "wall function" is used by the program when a specified wall is called for. 
Most walls would be input as specified walls with all velocities set to zero. A moving 
wall could be modeled by specifying non-zero velocities in the plane of the wall. A sym- 
metric boundary is a plane having no gradients across it. A periodic boundary is a 
boundary shared by both ends of the computational domain to simulate repeating geometries 
and flow patterns. Symmetric boundaries and periodic boundaries can be specified only for 
the outermost boundaries of the computational domain and must be specified identically for 
every t block bounding the symmetric or periodic plane. 

Section 3.4 gives a step by step procedure for specifying geometry input and the user is 
strongly advised to read it carefully before attempting to model his or her case. 

The following is a description of the card image input file. Note that any field which 
has a default value may be left blank. Caution should be observed when using default val- 
ues# however# because the default value may not be correct for the problem being solved. 
Also note that although English units are shown# 51 units may be substituted by setting 
IUNIT (card 7) to the appropriate value. 
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The format of each Input item is identified by type. 


TYPE 

A 

F 


I 


L 


earn s 


P5SCPIPTIQN 


0333:3333: 


ALPHANUMERIC ; Any keyboard characters are 
specified within the given field. 




b:333:::33 


FLOATING POINT : Decimal fractions including 
decimal point are specified anywhere within the 
given field. Positive values are assumed 
unless the value is preceded by a minus sign. 


Scientific notation may be substituted by 
specifying a decimal mantissa (as above) and a 
right adjusted base 10 multiplier preceded by 
a symbol E. 

INTEGER : Right adjusted whole numbers (no 

decimal point) are specified within the given 
f i eld. 

LOGICAL : T— true# F— false; right adjusted in the 

field. 


Blanks in an F or I field are interpreted as zeros# blanks in an A field are interpreted 
as blank characters# and a blank L field is interpreted as false. 


Card Type: 1 

Title: Title Card 

General Description: 


Name 

Type 

Columns 

Default 

Units 

Descri pt i on 

TITLE 

A 

1-80 



Descriptive title for case which will appear on 
printout and all plots. 
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Card Typo: 2 

Title: Origin and Countars 

General Description: 


Name 

Type 

Columns 

Default 

Uni ts 

Descri pt i on 

NXBLK 

I 

1-5 


Note 

1 

Number of blocks used in the block breakup 
X-di recti on . 

i n 

the 

NYBLK 

I 

6-10 


Note 

1 

Number of blocks used in the block breakup 
Y-di recti on. 

i n 

the 

NZBLK 

I 

11-15 


Note 

1 

Number of blocks used in the block breakup 
Z-di recti on . 

i n 

the 

NFBLK 

I 

16-20 


Note 

1 

Number of fluid block description cards 
type 4). 

( card 

NBC 

I 

21-25 


Note 

1 

Number of boundary condition description 
(card type 5A and 5B). 

cards 

IPAPER 

I 

26-30 

20 



Paper code for plots* See Note 2. 



XBEGIN 

F 

31-40 


ft 


Axial location of the south-west corner 
geometry. 

of 

the 

YBEGIN 

F 

41-50 


ft 


Radial location of the south-west corner 
geometry . 

of 

the 

ZBEGIN 

F 

51-60 


ft 


Tangential location of the front south-west cor- 
ner of the geometry. 

IROT 

I 

66-70 




Rotation Indicator 

IROT = 0 No rotation 

IROT = 1 Rotation in component vector 
form (card type 7C must be 
suppl i ed) 

IROT = 2 Rotation in magnitude and 
direction form (card type 
7D must be supplied) 




Notes : 

1. There is a limit on tho maximum number of these cards. See Appendix E. 

2. Appendix B has a listing of available paper codes. 


Card Type: 3A 

Title: West to East Block Break-Up 

General Description: 

This card repeated NXBLK times. 


Name 

Type 

Columns 

Default 

Uni ts 

Description 

XEND 

F 

1-10 


ft 

The X-di recti on coordinate of the east end of a 
plane of blocks. 

NX 

I 

11-15 



The number of cells the block is to be broken up 
into in the west to east direction. Note 1. 

EPSX 

■ 

16-20 

1.0 

Note 4 

The expansion factor to be used in this block. 
Expansion factor is the ratio of the width of a 
given cell to the width of its neighbor on the 
side closest to the origin. 
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ISSUED: 11/16/84 
REVISED: 03 / 04/86 

Card Typo: 3B 

Title: South to North Block Break-Up 

General Descri pt i on : 

Thi 5 card repeated NYBLK times. 


Name 

Type 

Columns 

Default 

Units 

Description 

YEND 

F 

1-10 


ft 

The Y-di recti on coordinate of the north end of a 
plane of blocks. 

NY 

I 

11-15 



The number of cells the block is to be broken up 
into in the south to north direction. Note 2. 

EPSY 

F 

0 
CM 

1 

*£> 

H 

1.0 

Note 4 

The expansion factor to be used in this block. 


Card Typo: 30 

Title: Front to Back Block Break-Up 

General Description: 

This card repeated NZBLK times. 


Name 

Type 

Columns 

Default 

Units 

Description 

ZEND 

F 

1-10 


ft 

The Z-di recti on coordinate of the back end of a 
plane of blocks. 

NZ 

I 

11-15 



The number of cells the block i s to be broken up 
into in the front to back direction. Note 3. 

EPSZ 

F 

16-20 

1.0 

Note 4 

The expansion factor to be used in this block. 


Notes: 

1* The sum of the NX’s must not exceed a maximum. See Appendix E. 

2. The sum of the NY’s must not exceed a maximum. 

3. The sum of the NZ f s must not exceed a maximum. 

4. Inaccuracies and instabilities can occur for values different than 1.0. Values greater 

than 1.25 or less than 0.8 should be avoided whenever possible. 
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ISSUED: 11/16/84 
REVISED: 03 / 04/86 

Card Typo: 4 

Title: Fluid Block Description 

General Description: 

This card repeated NFBLK times. 


Name 

Type 

Columns 

Default 

Units 

Descri pt i on 

MFBI 

I 

1-5 




Considered fluid block number in the X-direction. 

MFBJ 

I 

6-10 




Consider fluid block number in the Y-direction. 

MFBK 

I 

11-15 




Considered fluid block number in the Z-di recti on. 

MBCW 

I 

16-20 

Note 

i 


Integer name of boundary condition description 
cards used to describe the boundary conditions on 
the west boundary of this block (cards 5A and 
5B ) . 

MBCE 

I 

21-25 

Note 

1 


Integer name of boundary condition description 
cards used to describe the boundary conditions on 
the east boundary of this block (cards 5A and 
5B ) . 

MBC5 

I 

26-30 

Note 

1 


Integer name of boundary condition description 
cards used to describe the boundary conditions on 
the south boundary of this block (cards 5A and 
5B ) . 

MBCN 

I 

31-35 

Note 

1 


Integer name of boundary condition description 
cards used to describe the boundary conditions on 
the north boundary of this block (cards 5A and 
5B) . 

MBCF 

I 

t 

36-40 

Note 

1 


Integer name of boundary condition description 
cards used to describe the boundary conditions on 
the front boundary of this block (cards 5A and 
5B) 

MBCB 

I 

41-45 

Note 

1 

I 

Integer name of boundary condition description 
cards used to describe the boundary conditions on 
the back boundary of this block (cards 5A and 5B) 


Note 1: These integers can be anything between 1 and the maximum for NBC (Card Type 2) 

and are given in the first five columns of the Card Type 5A. Please see Appendix 
B> Sample Input/ for further clarification. 
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ISSUED: 11/16/84 
REVISED: 03 / 04/86 


Card Typo: 5A 

Title: Boundary Condition Description Card 

General Descri pt i on : 

There is one card 5A followed by one card 5B for each boundary condition. There are NBC 
card pairs to be input. 


Name 

Type 

Columns 

Default 

Units 

Descr i pt i on 

LBC 

I 

1-5 



Unique integer "name" assigned to this boundary 






condition card pair. 


IBC 

A 

8-10 



Boundary condition type: 







IBC = SYJ1 : axis of symmetry or 







plane of symmetry 

N 






IBC = UUL : unspecified wall 

o 






IBC = SUL : specified wall 

t 






IBC = UOP : unspecified opening 
IBC = SOP : specified opening 

e 






IBC = PER : periodic boundary 

1 

UBC 

F 

11-20 


ft/sec 

Velocity in X-direction at this boundary. 

VBC 

F 

21-30 


ft/sec 

Velocity in Y-direction at this boundary. 

UBC 

F 

31-40 


ft/sec 

Velocity in Z-di recti on at this boundary. 

PSBC 

F 

41-50 


lb/in 2 

Static pressure at this boundary. 


TBC 

F 

51-60 


°F 

Temperature at this boundary. Note 2. 

TEBC 

F 

61-70 



Turbulence energy coefficient at 
Note 3. 

this boundary. 

IUUNIT 

I 

71-75 



Swirl indicator at this boundary. 

Note 4. 






IUUNIT = 0 No conversion of UBC. UBC in ft/sec 






IUUNIT > 0 Solid body rotation. 







UBC input in rev/mi n 


EMISSW 

F 

76-80 



Emissivity of this boundary. 



Notes : 

1. The boundary condition types can also be input as right adjusted integers; 1 - SYM, 2 
=_UUL, 3 = SUL, 4 = UOP, 5 = SOP, 6 = PER. Also note that periodic and symmetry con- 
ditions can only be applied to the outer most boundaries and must be specified for 
every block bounding the symmetric or periodic plane. 

2. A value of -10,000. for temperature indicates to the program that the wall is adia- 
batic. 

3. The coefficient is the constant of proport i onal i ty in the equation: K = TEBC # U 2 

where TEBC = 3/2* (turbulence intensity) 2 

4. The units depend on the swirl indicator IWUNIT, 
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ISSUED: 11/16/84 
REVISED: 03/04/86 


Card Typo: 5B 

Title: Boundary Condition Description Card 

General Description: 

There is one card 5A followed by one card 5B for each boundary condition. There are NBC 
card pairs to be input. 


Name 

Type 

Columns 

Default 

Uni ts 

Descr i pt i on 

FBC(l) 

F 

1-10 



Mass fraction of species 1. Note 1. 

FBCC2) 

F 

11-20 



Mass fraction of species 2. Note 1 

FBCC3) 

F 

21-30 



Mass fraction of species 3. Note 1 

FBC(4) 

F 

31-40 



Mass fraction of species 4. Note 1. 

EDBC 

F 

71-80 


ft 2 / 

sec 3 

Energy disipation at this boundary. Note 2. 


Notes: 

1. A value of -10,000. for any mass fraction indicates to the program that there is no 
diffusion of the species across the boundary. 

2. The energy dissipation at all specified openings T S0P f must be specified. In the ab- 
sence of a better value, the following formula can be used EDBC = 3.0 X CTEBC) XX 1.5 

U 3 /L where TEBC was specified on card type 5A and L is the height of the passage. 
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ISSUED: 11/16/84 
REVISED: 03/04/86 

Card Type: 6 

Title: Case Indicators 

General Description: 


Name 

Type 

Columns 

Default 

Units 

Descri pt i on 

INDCOS 

I 

1-5 

2 


Problem type indicator 






INDCOS = 1 : Cartesian coordinates 
INDCOS = 2 : Cylindrical coordinates 

IPREF . 

I 

6-10 

1 


Number of cells in the X-di recti on of the refer- 
ence cell used for initial temperature and pres- 
sure guess. 

JPREF 

I 

11-15 

1 


Number of cells in the Y-di recti on of the refer- 
ence call used for initial temperature and pres- 
sure guess. 

KPREF 

I 

16-20 

1 


Number of cells in the 2-direction of the refer- 
ence cell used for initial temperature and pres- 
sure guess. 

IMON 

I 

21-25 

2 


Number of cells in the X-di recti on of cell to be 
monitored in the output file. 

JMON 

I 

26-30 

2 


Number of cells in the Y-direction of cell to be 
monitored in the output file. 

KMON 

I 

31-35 

2 


Number of cells in the Z-direction of cell to be 
monitored in the output file. 

DAMPW 

F 

41-50 

100. 

Note 1 

Swirl damping factor. 

DAMPR 

F 

51-60 

100. 

Note 1 

Radiation damping factor. 

TBMN 

F 

61-70 

1 . 

°R 

Minimum temperature limit for temperature predie- 
t i on . 

TBMX 

F 

71-80 

1.E6 

*R 

Maximum temperature limit for temperature predic- 
tion. 


Note 

1* These damping factors have been optimized to give optimum performance over a wide range 
of flow conditions and should be adjusted if a number of parametric studies are to bo 
conducted. 
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ISSUED: 1 1/16/8$ 
REVISED: 03/04/86 

Card Typo: 7 

Title: Case Indicators 

General Description: 


Name 

Type 

Columns 

Default 

Uni ts 

Descr i pt i on 

INLAM 

I 

1-5 

0 


Laminar flow indicator 

INLAM = 0 : laminar flow 
INLAM = 1 : turbulent flow 

IDENS 

I 

6-10 

0 


Density calculation indicator 

IDENS < 0 : Density calculations proceed every 

iteration starting with the -IDENSth 
i terat i on 

IDENS = 0 : Constant density case 

IDENS > 0 : Density calculations performed 
every IDENSth iteration 

INPRO 

l 

11-15 

F 

Note 1 


Fluid Property calculation switch 

INPRO = F : properties calculated every iter- 
ation 

INPRO = T : VISCOS (card 8) used for viscosity 
DEN5IT (card 8) used for density 

INDPRI 

I 

16-20 

1000 


All variable fields printed every INDPRI iter- 
ation. 

IUNIT 

I 

21-25 

1 


Input and output units indicator 
IUNIT = 1 : 51 units 

IUNIT = 2 : English units (lengths in feet) 
IUNIT = 3 : English units (lengths in inches) 

MAXIT 

I 

26-30 

200 


Maximum number of iterations allowed to converge 
before stopping. 

SORMAX 

F 

31-40 

. 005 


The maximum non di mensi onal i zed value each resi- 
dual source may have# to consider the case con- 
verged. 

I TEMP 

I 

41-45 

0 


Temperature calculation indicator 

ITEMP < 0 : Temperature calculations proceed 
every iteration starting with the 
-ITEMPth iteration 

ITEMP - 0 : Constant temperature case 
ITEMP > 0 : Temperature calculations performed 
every ITEMPth iteration 


Notes: 


1 . 


Fluid properties include eddy viscosity which has to bo calculated every iteration if 
turbulent flow is being calculated. Hence enter (INPR0=F) if INLAM>0. Also for vari- 
able density flows enter (INPR0=F). As a rule, if in doubt set (INPRQ=F). 
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Card Typo: 7 (continued) 


ISSUED: 1 1/16/84 
REVISED: 03/04/86 


Name 

Type 

Columns 

Default 

Uni ts 

Descr i pt i on 

KMBUST 

I 

46-50 



Combustion model indicator 






KMBUST = 0 : No combustion 
KMBUST = 1 : Eddy breakup model 

IRAD 

I 

51-55 

-50 


Radiation calculation indicator 






IRAD < 0 : Radiation calculations proceed every 
iteration starting with the -IRADth 
i terat i on 

IRAD > 0 : Radiation calculations performed 
every IRADth iteration 

EGAS 

F 

61-70 



Emmisivity of the gas. Overrides the emmisivity 

calculated by the program. 

DH 

F 

71-80 



Hydraulic diameter used for radiation 

calculati on . 


Card Type: 7C 

Title: Rotation in Component Vector Form 

General Description: 

This card describes the rotation of the control volume about the origin of the axes by 
giving the component vectors along the principal axes of the coordinate system. It is in- 
put when IR0T=1 (Card Type 2). 


Name 

Type 

Columns 

Default 

Units 

Descr i pti on 

OMEGAX 

F 

1-10 


Rad/sec 

Component 
X-axi s. 

of 

the 

rotat i on 

vector 

along 

the 

OMEGAY 

F 

11-20 


Rad/sec 

Component 
Y-axi s . 

of 

the 

rotat i on 

vector 

along 

the 

OMEGAZ 

F 

21-30 


Rad/sec 

Component 
Z-axi s . 

of 

the 

rotat i on 

vector 

along 

the 
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ISSUED: I 1/16/84 
REVISED: 03 / 04/86 


Card Type: 7D 

Titlo: Rotation in Magnitude and Direction 

General Description: 

This card describes the rotation of the control volume about the origin of the axes by 
giving the magnitude direction of the rotation vector relative to the principal axes. It 
is input only when IR0T=2. 


Name 

Type 

Columns 

Default 

Uni ts 

Descri pt i on 

0WGR0T 

F 

1-10 


Rad/sec 

Magnitude of the rotation of the control volume 
about the origin. 

THTROT 

F 

11-20 


Radi ans 

The angle between a projection of the roation 
vector on the XY-plane and the X-axis. 

PHIROT 

F 

21-30 


Radi ans 

The angle between the rotation vector and the 
Z-axi s. 


Card Type: 8 

Title: Fluid Properties 

General Description: 


Name 

Type 

Columns 

Default 

Units 

Description 

ALAMDA 

F 

1-10 

.005 


Length scale factor used to normalize energy dis- 
sipation residuals. 

VISCOS 

F 

11-20 


lbm/ 
f t-sec 

Fluid viscosity. 

DENSIT 

F 

ro 

H 

1 

O 


lbm/ft 3 

Density for constant density cases. 

FW 

F 

31-40 

29. 


Molecular weight of fluid for non-mixing cases. 

PRANDL 

F 

41-50 

.7 


Prandtl number for laminar flow cases. 

C3 

F 

51-60 

Note 1 

Constan 

of the turbulenco model. 

DECLIN 

F 

61-70 

0.0 

Deg 

Angle of declination measured from the vertical. 
Used only when INBUOY = T (card 9). 

CPDUM 

F 

71-80 

.25 

Btu/ 

( 1 bm 0 F ) 

Constant pressure specific heat of fluid. 
Note 2. 


Notes: 

1 . Defaults to the value of Cl on card type 12. 

2. Any negative value indicates to the program the enthalpy cards are to be read in. 
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ISSUED: 11/16/84 
REVISED: 03/04/86 

Card Typo: 9 

Title: Calculation Switches 

General Description: 

All switches turn calculation on if T is input, and off if F is input. 


Name 

Type 

Columns 

Default 

Uni ts 

Descr i pt i on 

INCALU 

L 

1-5 

F 


U-velocity calculation switch. 

INCALV 

L 

6-10 

F 


V-velocity calculation switch. 

INCALU 

L 

11-15 

F 


U-velocity calculation switch. 

INCALP 

L 

16-20 

F 

Note 1 

Pressure calculation switch. 

INCALK 

L 

21-25 

F 

Note 2 

Turbulence kinetic energy calculation switch. 

INCALD 

L 

26-30 

F 

Note 2 

Energy dissipation calculation switch. 

INCALT 

L 

31-35 

F 


Temperature calculation switch. 

INCALF 

L 

36-40 

F 


Mixing calculation switch. 

INVD 

L 

41-45 

F 

| 

Divergence of velocity calculation switch. 

INBUOY 

L 

46-50 

F 


Buoyancy calculation switch. 

INRAD 

L 

51-55 

F 


Radiation calculation switch. 

ISKEW 

I 

56-60 

F 

Note 3 

Difference scheme selector: 0-hydrid, 1-SUDS, 

2-BSUDS2 . 

DBSKEU 

I 

61-65 



DBSKEU>0 bounding variables printed. 


Notes: 

1- This switch in addition to calculating pressure, solves the continuity equation as well 
and must be turned on even if constant pressure or constant density calculations are being 
made. As a rule, if in doubt turn this switch on. 

2. Turbulence energy and dissipation must be calculated if turbulent flow is being com- 
puted. Hence these two switches must be turned on if INLAM>0 on card type 7. 

3. BSUDS2 and SUDS are differencing schemes that are more accurate than hybrid but these 
can become unstable for some flows. If instability occurs, starting with a converged hy- 
brid solution helps sometimes. But in some cases, it may not be possible to obtain a con- 
verged solution with BSUDS2 or SUDS. 
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ISSUED: 11/16/84 
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Card Type: 10 

Title: Under-relaxation Factors 

General Description: 

Unless there is a very unusual circumstance these variables should be allowed to assume 
their default values. 


Name 

Type 

Columns 

Default 

Uni ts 

Descri pt i on 

URFU 

F 

1-5 

.5 

Note 1 

U-velocity calculation under-relaxat i on factor. 

URFV 

F 

6-10 

.5 

Note 1 

V-velocity calculation undei — relaxation factor. 

URFI4 

F 

11-15 

.5 

Note 1 

W-velocity calculation under-relaxation factor. 

URFP 

F 

6-20 

1.0 

Note 2 

Pressure calculation under-relaxat i on factor. 

URFK 

F 

21-25 

.7 


Kinetic energy calculation under-relaxat i on fac- 
tor . 

URFE 

F 

26-30 

.7 


Energy dissipation calculation undei — relaxation 
factor . 

URFT 

F 

31-35 

1.0 


Temperature calculation under-relaxat i on factor. 

URFVIS 

F 

36-40 

.7 


Viscosity calculation under-relaxat i on factor. 

URFDEN 

F 

41-45 

.3 

Note 3 

Density calculation under-relaxat i on factor. 

URFF 

F 

46-50 

.7 


Mixing calculation under-relaxation factor. 

URFDIV 

F 

51-55 

.8 


Divergence calculation under-relaxat i on factor. 

URFR 

F 

56-60 

1.0 


Radiation calculation under-relaxation factor. 


Notes : 

1. It is not advisable to increase these under-relaxat i on factors beyond the default va- 
lue. Reducing these to 0.3 usually helps to converge cases which are unstable during 
the iteration process. 

2. The default value should be reduced to achieve convergence. However? a large penalty 
in convergence rate is incurred if this value is reduced by a large amount? below 0.8 
for example. 

3. For combustion cases and where large density gradients exist? this value may be re- 
duced to 0.1. 
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ISSUED: 11/16/84 
REVISED: 03 / 04/86 


Card Type: 11 

Title: Number of Sweeps 

General Description: 

Number of times each field variable is swept by the tridiagonal matrix solver. Unless 
there is a very unusual circumstance these variables should be allowed to assume their de- 
fault values. 


Name 

Type 

Columns 

Default 

Un i ts 

Descri pt i on 

NSUIPU 

I 

1-5 

3 


Number of U-velocity calculation sweeps. 

NSUIPV 

I 

6-10 

3 


Number of V-velocity calculation sweeps. 

NSWPW 

I 

11-15 

3 


Number of W-velocity calculation sweeps. 

NSWPP 

I 

16-20 

5 


Number of Pressure calculation sweeps. 

NSWPK 

I 

21-25 

3 


Number of Kinetic energy calculation sweeps. 

NSWPD 

I 

26-30 

3 


Number of Energy dissipation calculation sweeps. 

NSWPT 

I 

31-35 

3 


Number of Temperature calculation sweeps. 

NSWPF 

I 

46-50 

3 


Number of Mixing calculation sweeps. 

NSWPR 

I 

56-60 

10 


Number of Radiation calculation sweeps. 


Card Type: 12 

Title: Turbulent Air Properties 

General Description: 

Unless there is a very unusual circumstance these variables should be allowed to assume 
their default values. 


Name 

Type 

Columns 

Default 

Units 

Descri pt i on 

CMU 

F 

1-10 

.09 



CD 

F 

11-20 

1.0 



Cl 

F 

21-30 

1.44 



C2 

D 

31-40 

1.92 


Constants of the turbulence model. 

CAPPA 

B 

4* 

1 

U1 

o 

.4187 



ELOG 

n 

51-60 

9.793 



PRTE 

F 

61-70 

1. 



PRANDT 

F 

71-80 

.5 
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ISSUED: 11/16/84 
REVISED: 03 / 04/86 


Card Typo: 13A 

Title: Specie Names 

General Description: 


Name 

Typo 

Columns 

Default 

Uni ts 

Descri pt i on 

SPENAMC 1 ) 

A 

3-10 

f H2 f 


Name associated with specie 1. 

SPENAMC2) 

A 

13-20 

CNi 

O 


Name associated with specie 2. 

SPENAMC 3) 

A 

23-30 

1 C02 1 


Name associated with specie 3. 

SPENAMC4) 

A 

33-40 

f H20 1 


Name associated with specie 4. 


Card Type: 13B 

Title: Specie Switches 

General Description: 

Specie ”1" will be included in the mixture calculation if ICOMP(L) is T. It will be ex- 
cluded if ICOMP(L) = F. 


Name 

Type 

Columns 

Default 

Units 

Descri pt i on 

IC0MPC1) 

L 

1-10 

F 


Calculation switch for specie 1. 

IC0MPC2) 

L 

11-20 

F 


Calculation switch for specie 2. 

IC0MPC3) 

L 

21-30 

F 


Calculation switch for specie 3. 

I COMP (4) 

l 

31-40 

F 


Calculation switch for specie 4. 

MIXTYP 

I 

71-75 

1 


Mixture calculation type indicator 

MIXTYP = 1 : Assumes no specie fluctuations. 

Note 1. 


Note: 

1. For the present, only one mixture calculation is available so 1 must be specified. 
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ISSUED: 11/16/84 
REVISED: 03/04/86 

Card Typo: 13C 

Title: Specie Molecular Weights 

General Description: 

Specify the molecular weight for each specie selected on card type 13B. 


Name 

Type 

Columns 

Default 

Uni ts 

Descri pt i on 

XMOLWT(l) 

F 

1-10 

2.016 


Molecular weight of specie 1. 

XM0LWTC2) 

F 

11-20 

32. 


Molecular weight of specie 2. 

XM0LWTC3) 

F 

21-30 

44.01 


Molecular weight of specie 3. 

XM0LWTC4) 

F 

31-40 

18.016 


Molecular weight of specie 4. 


Card Type: 13D 

Title: Specie Laminar Flow Prandtl Number 

General Description: 

Specify the laminar flow Prandtl number for each specio selected on card 13B. 


Name 

Type 

Columns 

Default 

Units 

Description 

PRANMLC1) 

F 

1-10 

1 . 


Laminar flow Prandtl number of specie 1. 

PRANMLC2) 

F 

11-20 

1 . 


Laminar flow Prandtl number of specie 2. 

PRANMLC3) 

F 

21-30 

1 . 


Laminar flow Prandtl number of specie 3. 

PRANML (4) 

F 

31-40 

1 . 


Laminar flow Prandtl number of specie 4. 
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ISSUED: 11/16/84 
REVISED: 03 / 04/86 

Card Typo: 13E 

Title: Specie Turbulent Flow Prandtl Number 

General Description: 

Specify the turbulent flow Prandtl number for each specie selected on card 13B. 


Name 

Type 

Columns 

Default 

Units 

Descri pt i on 

PRANMTC1) 

F 

1-10 

0.5 


Turbulent flow Prandtl number of specie 1. 
Note 1. 

PRANMTC2) 

F 

11-20 

0.5 


Turbulent flow Prandtl number of specie 2. 

PRANMT (3 ) 

F 

21-30 

0.5 


Turbulent flow Prandtl number of specie 3. 

PRANMTC4) 

F 

31-40 

0.5 


Turbulent flow Prandtl number of specie 4. 


Note: 

1. Turbulent Prandtl numbers for all species should be equal to these values and should 
also be equal to PRANDT on card type 12. 

Card Type: 14 

Title: Combustion 

General Description: 

This card is input only if MIXTYP > 0 (card 13B). 


Name 

Type 

Columns 

Default 

Uni ts 

Descri pt i on 

IFUEL 

I 

1-5 

1 


The number of the fuel species. 

I0X 

I 

6-10 

2 


The number of oxidant species. 

IPR0D1 

B 

11-15 

3 


The number of the first product species. 

IPR0D2 

■ 

16-20 

4 


The number of the species to be the second prod- 
uct . 

STOIC 

F 

31-40 

7.937 


Stoichiometric oxidant/fuel ratio. The default 
value is for hydrogen and oxygen. 

FPR0D1 

F 

41-50 

1 . 


Fraction of first product in total products. 

A1MAG 

F 

51-60 

4. 


Combustion model constant. 

A2MAG 

F 

61-70 

2. 


Combustion model constant. 
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ISSUED: 11/16/84 
REVISED: 03/04/86 

Card Typo: 15 

Title: Mixture Fraction 

General Description: 

This card is input only if WIXTYP > 1 (card 13B). This card is used at present to define 
the value of the mixture fraction at the boundary. 


Name 

Type 

Columns 

Default 

Uni ts 

Descri pt i on 

IFUEIS 

1 

1-5 



Cell number x-di recti on for the fuel stream. 

JFUELS 

I 

6-10 



Cell number y-direction for the fuel stream. 

KFUEL5 

I 

11-15 



Cell number z-direction for the fuel stream. 

IOXS 

B 

16-20 



Cell number x-direction for the oxidant stream. 

J0X5 

B 

21-25 



Cell number y-direction for the oxidant stream. 

KOXS 

B 

26-30 



Cell number z-direction for the oxidant stream. 

PRFMl 

B 

31-40 



Laminar Prandtl number. 

PRFMT 

B 

41-50 



Turbulent Prandtl number. 

PRFLUC 

B 

51-60 



Prandtl number for fluctuation. 

CGI 

I 

61-70 

! 


Constant for mixture fraction fluctuation 






equation. 

CG2 

F 

71-80 



Constant for mixture fraction fluctuation 






equat i on . 
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3.3.2 Enthalpy Card Input 

The following is a description of the enthalpy card file. This file is used to set up a 
table of enthalpies vs. temperatures for each species involved in the case. These cards 
are required only if temperatures are being calculated and if CPDUM (card type 8) is nega- 
tive. The mixing calculation switch* INCALF (card type 9) must be turned on. 


Card Type: HI 

Title: Enthalpy Input - Species Names 

General Description: 


Name 

Type 

Columns 

Default 

Units 

Descri pt i on 

IUNTAB 

I 

1-5 

1 


Input units indicator for enthalpies 








IUNTAB = 1 : SI units 
IUNTAB = 2 : English units 




HNAMES(l) 

A 

11-18 



Name of 
plied. 

specie 1 for which enthalpy 

i s 

bei ng 

sup- 

HNAME5(2 ) 

A 

21-28 



Name of 
plied. 

specie 2 for which enthalpy 

i s 

being 

sup- 

HNAMES(3) 

A 

31-38 



Name of 
plied. 

specie 3 for which enthalpy 

i s 

bei ng 

sup- 

HNAMES(4) 

A 

41-48 



Name of 
plied. 

specie 4 for which enthalpy 

i s 

bei ng 

sup- 
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Card Type: H2 

Title: Enthalpy Input - Temperature and Enthalpies 

General Description: 

This card must be repeated as many times (up to 20) as are required to adequately describe 
the enthalpies of each specie as a function of temperature over the range of temperatures 
calculated for the case. 


Name 

Type 

Columns 

Default 

Units 

Description 

TTABLE 

F 

1-10 


°R 

Temperature for which enthalpies are being input. 

HTABLE(l) 

F 

11-20 


BTU/ 

(lb- 

mole) 

Enthalpy for species 1. 

HTABLEC25 

F 

21-30 


BTU/ 

(lb- 

mole) 

Enthalpy for species 2. 

HTABLEC3) 

F 

31-40 


BTU/ 

(lb- 

mole) 

Enthalpy for species 3. 

HTABLEC4) 

F 

41-50 


BTU/ 

(lb- 

mole) 

Enthalpy for species 4. 
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3.4 A Sample Input 

A step-by-step procedure for geometry Input is presented. For this purpose# a simple ge 
ometry is chosen (Figure 11), 


TURBULENT FLOW IN GEOFF'S DOGHOUSE GEOMETRY 




1.2 

0.8 


1.0 - 

- 1.2 



FIGURE 11 

Cut Auay Of Sample Geometry 
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It should be noted that TEACH models the fluid volume contained by the geometry shown in 
the previous figure. If the fluid volume is modeled with a coarse breakup, Figure 12 will 
result. It is clear from the modeling of this simple geometry that a preprocessor is re- 
quired to provide input to 3D-TEACH. In the absence of such a processor, a number of hand 
drawings and a lot of imagination will be needed to model even simple geometries in TEACH. 
From Figure 12 it is clear that five blocks are needed in the X-direction, four blocks are 
required in the y direction and three block are necessary in the Z-direction. 



FIGURE 12 

Block Breakup of Sample Geometry 
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Once the geometry shown in Figure 12 is visualized using the preprocessor or otherwise, 
the process of filling cards begins. 

Card No 1 is simple, we write TURBULENT FLOW IN GEOFFS DOGHOUSE. 


Card No 2 can only be filled in after cards up to 5 have been filled. 

Card 3A describes the cell break up. From Figure 12 we know the number of blocks in the 
X-direction, which gives us the number of cards of type 3A. Knowing the X-cordinate at 
which each of these blocks ends, also from Figure 12, XEND for all cards can be filled in. 
Now, the number of cells in each fluid block and the expansion factor associated with them 
have to be provided. This is dona by trail and error to give a smooth ovei — all break-up. 
For the present geometry, a uniform break-up is provided. 


XEND 


NX 


EPSX 


(3A) 


.4 

4 

1.0 

(3A) 


.5 

1 

1.0 

(3A) 


.8 

3 

1.0 

( 3A ) 


1.0 

2 

1.0 

( 3A ) 


1.2 

2 

1.0 

Cards 

3B 

and 3C can be 

filled similarly. 



YEND 

NY 

EPSY 

(3B) 


.5 

5 

1.0 

(3B) 


.7 

2 

1.0 

(3B) 


.8 

3 

1.0 

(3B) 


1.2 

4 

1.0 



ZEND 

NZ 

EPSZ 

(30 


.2 

2 

1.0 

(30 


.4 

2 

1.0 

(30 


.6 

2 

1.0 

Card ‘ 

typa 

4 should be 

filled in next; 1 

shows 

that it requires 

i the 

integer names 

vi ded 

i n 

card type 5. 

Hence 

card type 5 i 


however, inspection of card type 
of the boundary condition cards 
will be filled in first. 


4 description, 
which are pro- 


Figure 11 shows that there are only three different boundary conditions for this case, 
namely specified opening, wall, and unspecified openings. So three boundary condition 
card pairs (cards 5A and 5B, listed alternately) will be needed. The integer names (LBC) 
given to each card pair is arbitrary: 



LBC 

IBC 

UBC 

VBC 

WBC 

PSBC 

TBC 

TEBC 

(5A ) 

5 

SOP 

93. 

0 . 

0 . 

15. 

70 . 

.01 

#(5A) 

4 

UOP 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

( 5A ) 

3 

SWL 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 


FBC(l) 


FBC( 2 ) 

FBC( 3 ) 

FBC(4) 

EDBC 



(5B) 

0 


0 

0 

0 

24131 

. 0 . 


(5B) 

0 


0 

0 

0 

0 



( 5B ) 

0 


0 

0 

0 

0 




*Note: An unspecified opening can not be used in conjunction with an adjoining blockage. 
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For card 4, the number of fluid blocks that need to have boundary conditions applied has 
to be determined. This number can only be determined by filling in the card type 4 cards 
until the complete geometry is defined and then counting the number of cards filled. 
Filling these cards without the help of a pre-processor requires the geometry being mod- 
eled to be visualized in ones mind. Making several 2D sketches at different planes helps, 
but some trial and error is involved. In the absence of a pre-processor two tools are 
provided in the code. 1) The code checks these cards for consistent boundary conditions 
and flags the cards that are not consistent. This feature of the code is extremely useful 
in ensuring that the proper geometry is being specified. 2) Once the inconsistencies are 
corrected, the cell break up at different planes can be plotted to further check out the 
geometry. However these tools can only be used after the complete input file is set up. 
The recommendation at this point will be to fill these cards as best as one can, fill out 
other cards required for the INPUT, and proceed to run the case. Then using the two de- 
bugging tools described above, a complete and accurate speci f i cat i on of the desired 
geometry and mesh can be achieved. 

Once the iteration process described above is completed, it is found that the total number 
of cards required to specify the geometry is 39. For clarity, all the fluid blocks are 
shown in an exploded view in Figure 13. In this figure, block numbers are in parentheses 
and boundary condition numbers are circled. 
















TEACH 3D 


PAGE 56 


ISSUED: 11/16/84 
REVISED: Q3/04/86 


The "names" of the boundary condition description cards appear at all block boundaries. 
The boundary condition description cards for the six blocks annotated in Figure 13 are 
then completed: 


X-i ndex 

Y-i ndex 

Z- i ndex 

west 

east 

south 

north 

front 

back 

MFBI 

MFBJ 

MFBK 

MBCW 

MBCE 

MBCS 

MBCN 

MBCF 

MBCB 

1 

1 

1 

3 

3 

3 

0 

3 

0 

1 

1 

2 

5 

3 

3 

0 

0 

0 

1 

1 

3 

3 

3 

3 

0 

0 

3 

1 

2 

1 

3 

3 

0 

0 

3 

0 

1 

2 

2 

3 

3 

0 

0 

0 

0 

1 

2 

3 

3 

3 

0 

0 

0 

3 


For example* block (1*1*2) has boundary description cards number 5 describing its west 
boundary and 3 describing its east and south boundaries. A 0 on its north* front* and 
back boundaries indicates a boundary common with another fluid block. This process is re- 
peated until all 39 cards are filled in. 

This input is further explained in Figure 14 which illustrates the completed geometry in- 
put showing how boundary condition number 3 is applied to the west face of cell (1*1*1) 
from the breakup of Figure 12. 

Card Type 2 can now be filled. 

NXBLK NYBLK NZBLK NFBLK NBC IPAPER XBEGIN YBEGIN ZBEGIN 
5 4 3 39 3 

The paper code for Calcomp plots can be left blank to assume default value. The coordi- 
nates of lower south-west-f ront corner of the geometry can also be omitted* therby assum- 
ing default values* XBEGIN=0.0* YBEGIN=0.0, and ZBEGIN=0.0. 

Card Type 6: On this card* certain indicators and damping factors can be specified* but 
for the present problem* only the first four are needed. Since the geometry is cartesian* 
INDC0S = 1. It should be noted that IPREF* JPREF and KPREF should define the cel 1 * not 
block* which is located in the specified inflow boundary. Hence Card 6 can be filled in 
as follows: 

INCOS IPREF JPREF KPREF IMON JMON KMON DAMPW DAMPR TBMN TBMX 
112 3 

Card Type 7:^ This card allows one to run either a laminar or turbulent flow case and pro- 
vides the ability to use SI or English units. Maximum number of iterations* and the con- 
vergence level can also be specified. Ule would be running a turbulent case using English 
units for 200 iteration and the convergence level is 0.005. Hence card type 7 can be 
filled as follows. 

INLAfl IDENS INPRO INDPRI IUNIT MAXIT SORMAX ITEMP 
F F 2 200 0.005 
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SSOOG TURBULENT 
5 * 5 

0.4000 4 1. 

0.5000 1 1 

o.aooo 5 i 

1.0000 2 1 

1.2000 2 1 

0.5000 5 I 

0.7000 2 1 

o.aooo i i 

1.2000 4 1 

0.2000 2 1 

0.4000 2 1 

0.6000 2 I 

i i nr 


IN GEOFS DOGHOUSE 

5 10 0.0 


<33 £D 






l STM 

3) suV 
^ 0.0 

4 UOP 

0.0 

6 PER 

0.0 

5 SOP 

0.0 


f 0.0 

0.0 

0.0 

0.0 

o 

o 

0.0 

0 0.0 

0.0 

0.0 

0.0 




0. 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 1.00 

0.0 

0.0 

0.0 




0. 

0.0 

0.0 

0.0 

0.0 

70.00 

0.0 

0 1.00 

0.0 

0.0 

0.0 




0. 

0.0 

0.0 

0.0 

o 

o 

0.0 

0.0 

0 1.00 

0.0 

0.0 

0.0 




0. 

93.00 

0.0 

0.0 

15.00 

70.00 

0.01 

0 1.00 

0.0 

0.0 

0.0 




24131 . 


Figure 1^ 
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INDPRI* which specifies the number of iterations after which variable fields have to be 
printed* can be left to its default value. IDENS and ITEMP are not required for the pres- 
ent problem. 

Card Type 8: On this card* fluid properties are specified. For the present case* viscos- 
ity* density* molecular weight* and Prandtl number need to be specified. Other variables 
assume default values. 

ALAMDA VISCOS DENSIT FW PRANDL C3 DECLN CPDUM 
1.11E-5 0.0805 29 0.7 

Card Type 9: This card is used to specify the variable fields that need to be calculated. 
Since we are solving a case with turbulent flow* turbulence kinetic energy and its dissi- 
pation rate have to be calculated in addition to the three momentum equations and the con- 
tinuity equation. All other variables can be left to assume default values. 

INCALU INCALV INCALW INCALP INCALK INCALD 
T T T T T T 

Cord Type 10: This card specifies the under-relaxat i on factors. It was found by trail 
and error that the following under-relaxat i on factors worked best for the present config- 
urat i on. 

URFU URFV URFW URFP URFK URFE 
0.3 0.3 0.3 1.0 0.5 0.5 

This completes the number of INPUT cards that have to be filled. All the remaining cards 
can be left to assume default values and can be left blank. 


Figure 15 shows the complete INPUT file. 
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3.5 A Sample Output 


This section provides and explains a sample output. This output is a result of a TEACH 
run from the input given in the proceeding section. 





xi sysMMtrie 



Mracter i sties 


curist ically 


VERSION llO 
( 05 / 05 / 63 ) 




>j 



0.0 

0.0 

0.0 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 






9.0 


0.0 


0 


0 


35006 TURBULENT FLOW IN 6E0F3 
5 4 3 39 5 10 

0.4000 4 1.00 

0.5000 1 1.00 

0.8000 3 1.00 

1.0000 2 1.00 

1.2000 2 1.00 

0.5000 5 1.00 

0.7000 2 1.00 

0.0000 l 1.00 

1.2000 4 1.00 

0.2000 2 1.00 

0.4000 2 1.00 

0.4000 2 1.00 


DOGHOUSE 

0.0 


1 

1 


3 

3 3 

0 

3 

9 




4 

1 

2 


3 

3 0 

0 

3 

0 




4 

1 

3 


3 

3 0 

0 

3 

0 




4 

1 

4 


3 

0 0 

3 

3 

0 




4 

2 

4 


0 

0 3 

3 

3 

0 




4 

3 

4 


0 

0 3 

3 

3 

0 




4 

4 

4 


0 

3 0 

3 

3 

0 




4 

4 

3 


3 

3 0 

0 

3 

0 




4 

4 

2 


0 

3 0 

0 

3 

0 




4 

3 

2 


3 

0 0 

3 

3 

0 




4 

3 

1 


3 

0 3 

0 

3 

0 




4 

4 

1 


0 

0 3 

0 

3 

0 




4 

5 

1 


0 

4 3 

3 

3 

0 




4 

1 

1 


5 

3 3 

0 

0 

0 




4 

1 

2 


3 

3 0 

0 

0 

0 




4 

1 

3 


3 

3 0 

0 

0 

0 




4 

1 

4 


3 

0 0 

3 

Q 

0 




4 

2 

4 


0 

0 3 

3 

0 

0 




4 

3 

4 


0 

0 3 

3 

0 

9 




4 

4 

4 


0 

3 0 

3 

0 

0 




4 

4 

3- 


3 

3 0 

0 

0 

0 




4 

4 

2 


0 

3 0 

0 

0 

0 




4 

3 

2 


3 

0 0 

3 

0 

9 




4 

3 

1 


3 

0 3 

0 

0 

0 




4 

4 

1 


0 

0 3 

0 

0 

0 




4 

5 

1 

2 

0 

4 3 

3 

0 

0 




4 

1 

l 


3 

3 3 

0 

0 

3 




4 

1 

2 


3 

3 0 

0 

0 

3 




4 

1 

3 


3 

3 0 

0 

0 

3 




4 

1 

4 


3 

0 0 

3 

0 

3 




4 

2 

4 


0 

0 3 

3 

0 

3 




4 

3 

4 


0 

0 3 

3 

0 

3 




4 

4 

4 


0 

3 0 

3 

0 

3 




4 

4 

3 


3 

3 0 

0 

0 

3 




4 

4 

2 


0 

3 0 

0 

0 

3 




4 

3 

2 


3 

0 0 

3 

0 

3 




4 

3 

1 


3 

0 3 

0 

0 

3 




4 

4 

1 


0 

0 3 

0 

0 

3 




4 

5 

l 


0 

4 3 

3 

0 

3 




4 

1 

SYH 


0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

0 0.0 

5A 


0.0 


0.0 

0.0 


0.0 




0. 

58 

3 

SUL 


0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

0 1.00 

5A 


0.0 


0.0 

0.0 


0.0 




0. 

58 

4 

UOP 


0.0 

0.0 


0.0 

0.0 

70.00 

0.0 

0 1.00 

5A 


0.0 


0.0 

0.0 


0.0 




0. 

5B 

4 

PER 


0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

0 1.00 

5A 


0.0 


0.0 

0.0 


0.0 




0. 

58 

5 

SOP 

93.00 

0.0 


0.0 

15.00 

70.00 

0.01 

0 1.00 

5A 


0.0 


0.0 

0.0 


0.0 
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5B 

1 

1 

2 

3 

25 8 

0 


100. 

100. 

1. 

1000000. 

6 

F 

0 

F 

1000 

2 200 

0. 

0050 

0 0 

0 

-1.00000 

0.0 

7 

0. 

0050 . 

1I10E-04 

,80502-01 

29.00 

0.700 

1.44000 

0.0 

0.25000 

8 

T 

T 

T 

T 

T T 

F 

F 

F F 

F 



9 

000 

1.3000, 

3001 

.0000 

.3000.3001, 

0000.7000. 

3000.7000.8001.000 



10 

3 

3 

3 

5 

3 3 

3 

3 

10 




11 

0. 

0900 

1. 

0000 

1.4400 

l 

9200 

0.4167 

9.7930 

1.0000 

0.9000 

12 












13A 


F 


F 

F 


F 




1 

138 


2.01400 
1.00000 
0.90000 
1 2 


32.00000 
1.00000 
0.90000 
3 4 


44.00999 

1.00000 

0.90000 


18. 01401 
1.00000 
0.90000 
7.9370 


1.0000 


4.0000 


2.0000 


2 

3A 

3A 

3A 

3A 

3A 

38 

38 

38 

38 

3C 

3C 

3C 


13C 

130 

13E 

14 
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The next few pages of the output give information on the type of case being run. An exam- 
ple of this is given in Figure 18. 

Turbulent . 

Constant density. 

Viscosity is being calculated every iteration. 

All variable fields are being printed every 1000 iterations. 

All input in English units (lengths in 

All calculations sill stop after 200 iterations or when all residual sources fell below .50000E-02. 

No cowbustion. 

Ceeisivity of the gas *-1.0000 over ides the eesisivity calculation by the prograe. 

Rectangular coordinates. 

Calculations turned on* 

U- Velocity 
V-Veloci ty 
H- Velocity 
Pressure 
Kinetic Energy 
Energy Disipation 

THOMAS” ALGORITHM IS USED TO SOLVE MATRIX EQUATIONS. 

41 ERROR* IPREF* 1» JPREF* 2, KPREF* 3, IS NOT A SPECIFIED FLUID CELL. 

ERROR: IPREF* If JPREF* 2. KPREF* At HILL BE USED. 

At x ■ .0 free Y * .0 to Y ■ .5000 and free Z * .2000 to Z * .4000 

Mass floe in .74065 (LBM/SCO 

Inlat fluid velocity 93.000 (FT/S) 

Inlet fluid temperature 70.000 (DEG F) 

Reynolds number .26978E+06 

Fluid density .00500E-01 (LB/FT*P3) 

AT X * 0.0 FROM T* 0.0 TO Y * 0.50000 AND FROM Z ■ 0.20000 TO Z * 0.40000 

THE SUM OF THE FLOWS FOR HEST HALL IS 0. 74065 

Prandtl number ■ 7.000E-01 


Total ness floe in .74065 (UH/SEC) 

Cards have been raad in 

At X * .0 froe Y * .0 te Y * .5000 and free Z * .2000 to Z * .4000 

Mass floe in .74065 ( LBM/3EC) 

Inlat fluid velocity 93.000 (FT/S) 

Inlat fluid temperature 70.000 (DEG F) 

Reynolds number .26978E»06 

Fluid density .00500E-01 (LB/FT**3) 

AT X * 0.0 FROM Y* 0.0 TO Y * 0.50000 AND FROM Z * 0.20000 TO Z * 0.40000 

THE 5UH OF THE FLOWS FOR HE3T HALL IS 0.74065 

Prandtl nwbar ■ 7.000E-0! 


Total ease floe in .74065 ( LBM/SEC) 

Figure 18 

useri called. Case Information 

Note the error# *# in the center of the page. TEACH checks your boundary condition and if 
possible corrects small errors. All errors of this type will be noted here in the output. 
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Figure 19 below gives the next section of the output. Here the residuals of each field 
variable for each iteration are printed out. 


ITER I— ABSOLUTE RESIDUAL SOURCE 


HO. UMOM 

VHOW 

MASS 

TKIN TEMP 

1 .500 

.335E-05 

2.00 

. 102E-01 .0 

2 1.11 

.549 

2.30 

.2132-01 .0 

3 1.42 

2.22 

2.57 

.2452-01 .0 

4 2.14 

1.95 

1.60 

.239E-01 .0 

5 1.74 

2.03 

.875 

.2802-01 .0 

4 1.57 

1.99 

.527 

.3402-01 .0 

7 1.34 

1.75 

.513 

.3892-01 .0 

8 1.19 

1.55 

.528 

.4032-01 .0 

9 .997 

1.30 

.608 

.3912-01 .0 

10 .820 

1.02 

.595 

.3572-01 .0 

11 .750 

.926 

.557 

. 323E-01 .0 

12 .843 

1.00 

.499 

.2912-01 .0 

13 .923 

1.14 

.474 

.2702-01 .0 

14 .947 

1.18 

.324 

.2522-01 .0 

15 .901 

1.13 

.257 

.2372-01 .0 

14 .793 

1.01 

.273 

.2232-01 .0 

17 .644 

,854 

.334 

.2102-01 .0 

18 .500 

.659 

.314 

.1942-01 .0 

19 .392 

.502 

.272 

.1772-01 .0 

20 .402 

.471 

.230 

.1622-01 .0 

21 .447 

.498 

.192 

.1512-01 .0 

22 .470 

.510 

.153 

.1432-01 .0 

23 .457 

.495 

.127 

.1452-01 .0 

24 .406 

.439 

.116 

.1482-01 .0 

25 .328 

.349 

.115 

.1482-01 ,0 

26 .242 

.299 

.113 

.1442-01 .0 

27 .180 

.237 

.105 

.1341-01 .0 

28 .162 

.218 

.9O2E-01 

.1252-01 .0 

29 .170 

.228 

.712E-01 

.1122-01 .0 

30 .178 

.238 

.525E-01 

.9972-02 .0 

31 .176 

.240 

.3952-01 

.6942-02 .0 

32 .160 

.226 

. 396E-01 

.8172-02 .0 

33 .135 

.200 

.444E-01 

.7522-02 .0 

34 .106 

.166 

-471E-01 

.6922-02 .0 

35 .835E-01 

.134 

.4411-01 

.6352-02 .0 


34 .7872-01 .113 .4211-01 .5832-02 .0 

37 .841E-01 .109 .3562-01 .54IE-0Z .0 

38 .9112-01 .110 .2852-01 .5092-02 .0 

39 . 916E-01 .111 .2222-01 .4632-02 .0 

40 .849E-01 .105 .2062-01 .4642-02 .0 

41 . 737E-01 .9452-01 .2162-01 .450E-02 .0 

42 .603E-01 .SUE-01 .227E-01 .4351-02 .0 

43 .490E-01 .675E-01 .21SE-01 .4162-02 .0 

44 .424E-01 .595E-01 .193E-01 .394E-02 .0 

45 .4062-01 .5S9E-01 .164E-01 .369E-02 .0 

44 .416E-01 .594E-01 .134E-01 .344E-02 .0 

47 .405E-01 . 5SSE-01 .114E-01 .323E-02 .0 

48 ♦ 374E-01 .550E-01 .107E-01 .30SE-02 .0 

49 .330E-01 . 494E-01 .U0E-01 .293E-02 .0 

50 . 274E-01 .422E-01 .1132-01 .276E-02 .0 


SUMS 1 

DISP FLUC FRAC DEWS WHOM TR4 
.1772-03 .0 .0 .0 . 154E-04 .0 

. 113E-03 .0 .0 .0 .792 .0 

, 940E-04 .0 .0 ,0 .601 .0 

.835E-04 .0 .0 .0 .418 .0 

. 767E-04 .0 .0 .0 .345 .0 

. 722E-04 .0 .0 .0 .278 .0 

.6952-04 .0 .0 .0 .259 .0 

.678E-04 .0 .0 .0 .243 .0 

.664E-04 .0 .0 .0 .234 .9 

.648E-04 .0 .0 .0 .231 .0 

.630E-04 .0 .0 .0 .218 .0 

.610E-04 .0 .0 .0 .206 .0 

.601E-04 .0 .0 .0 .194 .0 

.588E-04 .0 .0 .0 .182 .0 

.5742-04 .0 .0 .0 .171 .0 

.556E-04 .0 .0 .0 .144 .0 

.533E-04 .0 .0 .0 .158 .0 

.508E-04 .0 .0 .0 .152 .0 

.483E-04 .0 .0 .0 .147 .0 

.468E-04 .0 .0 .0 .142 .0 

.462E-04 .0 .0 .0 .138 .0 

.459E-04 .0 .9 .0 .133 .0 

.457E-04 .9 .0 .9 .127 .9 

.457E-04 .9 .9 .0 .120 .0 

♦456E-04 .0 .0 .0 .112 .0 

. 449E-04 .0 .0 .0 .104 .0 

.437E-04 .0 .0 .0 .952E-01 .0 

.418E-04 .0 .0 .0 .849E-01 .9 

. 394E-04 .9 .9 .0 .792E-01 .9 

.3702-04 .0 .0 .9 .715E-01 .0 

.3442-04 .0 .0 ,9 .64IE-01 .9 

.317E-04 .9 .0 .0 .571E-01 .0 

» 289E-04 ,0 .0 .0 .510E-01 .0 

. 2622-04 .0 .0 .0 -455E-01 .0 

. 238E-04 .0 .0 .0 .4UE-01 .0 

.218E-04 .0 .0 .0 . 376E-01 .9 

. 202E-04 .0 .0 .0 . 348E-01 .9 

. 190E-04 .0 .0 .9 . 325E-01 .0 

. 182E-04 .9 .0 .0 . 303E-01 .0 

♦ 175E-04 .0 .0 .0 .284E-01 .0 

• 172E-04 .0 .0 .0 .265E-01 .0 

. 168E-04 .0 .0 .9 . 249E-01 .0 

. 164E-04 .0 .0 .9 .234E-01 .0 

.159E-04 .0 .0 ,9 . 221E-01 .0 

. 153E-04 .0 .0 .0 .209E-01 .0 

.146E-04 .9 .0 .0 . 198E-01 .0 

. 1 38E-04 .0 .0 .0 . 188E-01 .9 

. 130E-04 .0 .0 .0 . 179E-01 .0 

. 122E-04 .0 .0 .0 . 1 74E-01 .0 

. 1 14E-04 .0 .0 .0 . 147E-01 .0 


Figure 19 
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ITER I 

MO. UMOM 

VM0H 

-ABSOLUTE 

MASS 

RESIDUAL SOURCE 
TKIH TEMP 

SUMS — — — — 

DISP FLUC 

— I 

FRAC 

DENS 

HMOM 

TR4 

51 .2192-01 

.3482-01 

. 108E-01 

.2432-02 .0 

. 107E-04 .0 

.0 

.0 

. 162E-01 

.0 

52 . 188E-01 

.3082-01 

.9922-02 

. 248E-02 .0 

.1012-04 .0 

.0 

.0 

. 157E-01 

.0 

53 . 178E-01 

.305E-01 

.660E-0Z 

. 234E-02 .0 

.9432-05 .0 

.0 

.0 

.1522-01 

.0 

54 . 179E-01 

.315E-01 

.719E-02 

.2222-02 .0 

.9202-05 .0 

.0 

.0 

. 148E-01 

.0 

55 .180E-01 

. 318E-01 

.591E-02 

.21 IE-02 .0 

.0022-05 .0 

.0 

.0 

. 143E-01 

.0 

54 . 169E-01 

.309E-01 

.5492-02 

.203E-02 .0 

.0462-05 .0 

.0 

.0 

.137E-01 

.0 

57 . 151E-01 

.2862-01 

.5572-02 

. 195E-02 .0 

.0132-05 .0 

.0 

.0 

. 132E-01 

.0 

5* . 133E-0I 

. 255E-01 

.570E-02 

. 187E-02 .0 

.7792-05 ,0 

.0 

.0 

.1282-01 

.0 

59 . 1232-01 

. 2232-0] 

.557E-02 

. 176E-02 .0 

.7432-05 .0 

.0 

.0 

.1222-01 

.0 

40 . 122E-01 

• 198E-01 

.527E-02 

. 149E-02 .0 

.7072-05 .0 

.0 

.0 

.1172-01 

.0 

41 .1232-01 

. 187E-01 

.4632-02 

. 160E-02 .0 

.6692-05 .0 

.0 

.0 

.1122-01 

.0 

42 . 125E-01 

. 102E-01 

■429E-02 

, I50E-02 .0 

.630E-05 .0 

.0 

.0 

. 107E-OI 

.0 

43 .123E-01 

. 173E-01 

.379E-02 

. I42E-02 .0 

.5912-05 .0 

.0 

.0 

.1012-01 

.0 

44 . 117E-01 

. 16QE-Q1 

. 350E-02 

. 133E-02 .0 

.553E-05 .0 

.0 

.0 

.9632-02 

.0 

45 .1092-01 

. I43E-01 

. 341E-02 

. 125E-02 .0 

.5172-05 .0 

.0 

.0 

.914E-02 

.0 

44 .9642-02 

. 124E-01 

. 330E-02 

* HOE-02 .0 

.485E-05 .0 

.0 

.0 

.8452-02 

.0 

47 .8812-02 

. 1092-01 

. 308E-02 

.1112-02 .0 

.455E-05 .0 

.0 

.0 

.814E-02 

.0 

40 .706E-02 

. 101E-01 

. 280E-02 

.1052-02 .0 

.4292-05 .0 

.0 

.0 

. 764E-0Z 

.0 

49 .7232-02 

. 9882-02 

. 250E-02 

.100E-02 .0 

.4062-05 .0 

.0 

.0 

. 715E-02 

.0 

70 .470E-02 

.975E-02 

.2172-02 

.9412-03 .0 

.3072-05 .0 

.0 

.0 

.6682-02 

.0 

71 .445E-02 

.9562-02 

. 190E-02 

.9271-03 .0 

.3712-05 .0 

.0 

.0 

.626E-02 

.0 

72 .6072-02 

.9282-02 

.177E-02 

.0982-03 .0 

.3572-05 .0 

.0 

.0 

.588E-02 

.0 

73 .563E-02 

.8902-02 

.176E-02 

.0702-03 .0 

. 344 E- 05 .0 

.0 

.0 

.553E-02 

.0 

74 .521E-02 

.052E-02 

. 1762-02 

.042E-03 .0 

.3322-05 .0 

.0 

.0 

.524E-0Z 

.0 

75 .4902-02 

.8162-02 

. 174E-02 

.6142-03 .0 

.3222-05 .0 

.0 

.0 

•498E-02 

.0 

74 .473E-02 

.791E-02 

. 170E-02 

. 785E-03 .0 

.3122-05 .0 

.0 

.0 

.4762-02 

.0 

77 .4432-02 

. 700E-02 

, 163E-02 

. 754E-03 .0 

.3022-05 .0 

.0 

.0 

.4572-02 

.0 

70 .454E-02 

.776E-02 

. 153E-02 

.7282-03 .0 

.2932-05 .0 

.0 

.0 

.441E-02 

.0 

79 .440E-02 

.7692-02 

. 144E-02 

.7012-03 .0 

.2052-05 .0 

.0 

.0 

.4272-02 

.0 

00 .419E-02 

.7552-02 

. 145E-0Z 

.477C-03 .0 

.2772-05 .0 

.0 

.0 

.4132-02 

.0 

01 . 394E-02 

. 732E-02 

» 146E-02 

.4562-03 .0 

.2702-05 .0 

.0 


.400E-02 

.0 

02 . 347E-02 

. 704E-02 

. 147E-02 

.6372-03 .0 

.2642-05 .0 

.0 

.0 

. 388E-02 

.0 

03 . 342E-02 

.476E-02 

» 1452-02 

.6212-03 .0 

.2502-05 .0 

.0 

.0 

.375E-02 

.0 

04 .3222-02 

.6582-02 

. 143E-02 

,6072-03 .0 

.2532-05 .0 

.0 

.0 

.3442-02 


05 .3062-02 

.642E-02 

. J39C-02 

.5952-03 .0 

.2492-05 .0 

.0 

.0 

.353E-02 

.0 

04 . 2 94 E- 02 

-624E-02 

.1342-02 

.5042-03 .0 

. 245E-05 .0 

.0 

.0 

. 343E-02 

.0 

07 .282E-02 

.6042-02 

. 1292-02 

.5742-03 .0 

.2412-05 .0 

.0 

.0 

.332E-02 

.0 

00 .2692-02 

•570E-02 

♦ 1232-02 

.565E-03 .0 

.2382-05 .0 

.0 

.0 

.3202-02 

.0 

09 .2542-02 

.5492-02 

. 122E-02 

,5552-03 .0 

.2342-05 .0 

.0 

.0 

.3082-02 

.0 

90 .2392-02 

.516E-02 

. 120E-02 

.5442-03 .0 

.2302-05 .0 

.0 

.0 

.2962-02 

.0 

91 . 2242-02 

.4O4E-02 

. 1172-02 

.5332-03 .0 

.2272-05 .0 

.0 

.0 

.284E-02 

.0 


Figure 19 (continued) 

The rest of the output file consists of a listing of all field variables on all planes. 
The entire U-Velocity is given in Figure 20. The orginal output would also contain the 
following field variables: V-Velocity, W-Velocity, Static Pressure* Turbulence Energy, En- 
ergy Dissipation Rate, Viscosity, Density, East-West Peclet Number, North-South Peclet 
Number, Front-Back Peclet Number, Wall Reynolds Number, West and East Wall Shear Stress, 
South and North Wall Stress, and Front and Back Wall Stress. 
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U VELOCITY 


AT Zt 1 1*- .50000E-01 



J 

T 

2 

3 

4 

5 

6 

7 

6 

9 

14 

1.25000 

0.0 

0.0 

0.9 

0.0 

0.0 

0.0 

0.0 

0.0 

13 

1.15000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

12 

1.05000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

il 

0.95000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

10 

0.45000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

9.9 

0.0 

9 

0.75000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 

0.45000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

7 

0.55000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 

0.45000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

5 

0.35000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

4 

0.25000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3 

0.15000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

2 

0.05000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 

-0.05000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

X* 


0.0 

0.10000 

0.20000 

0.30000 

0.40000 

0.50000 

0.40000 

0.70000 

J 

Y 

10 

11 

12 

13 

14 




14 

1.25000 

0.0 

0.0 

0.0 

0.0 

0.0 




13 

1.15000 

0.0 

0.0 

0.0 

0.0 

0.0 




12 

1.05000 

0.0 

0.0 

0.0 

o.o 

0.0 




11 

0.95000 

0.0 

0.0 

0.0 

0.0 

0.0 




10 

0.65000 

0.0 

0.0 

0.0 

0.0 

0.0 




9 

0.75000 

0.0 

0.0 

0.0 

0.0 

o.o 




a 

0.450Q0 

0.0 

0.0 

0.0 

0.0 

0.0 




7 

0.55000 

0.0 

0.0 

0.0 

0.0 

0.0 




4 

0.45000 

0.0 

0.0 

0.0 

0.0 

0.0 




3 

0.35000 

0.0 

0.0 

0.0 

0.0 

0.0 




4 

0.25000 

0.0 

0.0 

0.0 

0.0 

0.0 




3 

0.15000 

0.0 

0.0 

0.0 

0.0 

0.0 




2 

0.05000 

0.0 

0.0 

0.0 

0.0 

0.0 




1 

-0.05000 

0.0 

0.0 

0.0 

0.0 

0.0 




X* 


0.60000 

0.90000 

1.00000 

1.10009 

1.20000 




»-*- 



U VELOCITY 



AT Zl 21* 

•50000E-01 



J 

Y 

2 

3 

4 

5 

4 

7 

6 

9 

14 

1.25000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

13 

1.15000 

0.0 

4.33500E+00 

1.047771+01 

1 .60506E+01 

3. 167252+01 

3.65556E+01 

3.m04E+0I 

3.70363E+01 

12 

1.05000 

0.0 

-1 . 17944E-01 

2.67547E+00 

1.24496E+01 

2 . 92490E+01 

3.53063E+01 

3.42420E+01 

3. 39066E+01 

11 

0.95000 

0.0 

-5.59304E+00 

-5.39041E+00 

7. 14951E+00 

3.04696E+01 

3.39939E+01 

3.50799E+01 

3.60562E+01 

10 

0.65000 

0.0 

-1 . 13273E+0I 

-1 . 35052E+01 

-4.76053E-02 

3.28602E+01 

2.42321E+01 

2 . 55091E+01 

3.44220E+91 

9 

0.75000 

0.0 

-1.69191E+01 

-2.202862+01 

-1.41953E+01 

0.0 

0.0 

0.0 

0.0 

a 

0.45000 

0.9 

-2. 11613E+01 

-2. 78947E+01 

-1 . 99235E+01 

0.0 

0.0 

1.O20I6E+Q1 

1 . 04476E+01 

7 

0.55000 

0.0 

-2.42204E+01 

-3. 17203E+01 

-2 . 27362E+01 

0.0 

0.0 

6. 14431E+00 

1 .02752E+0 1 

4 

0.45000 

0.0 

-2.63003E+01 

-3.37121E+01 

-2.36429E+01 

0.0 

0.0 

4.63557E+00 

1.23152E+01 

5 

0.35000 

0.0 

-2.639892+01 

-3.41044E+01 

-2. 36244E+0I 

0.0 

0.0 

4.64790E+00 

1 .06192E+01 

4 

0.25000 

0.0 

-3. 09010E+01 

-3.44254E+01 

-2.41 161E+01 

0.0 

0.0 

5. 73944E-01 

1 . 90141E+00 

3 

0.15000 

0.0 

-3.26551E+01 

-3.53937E+01 

-2.44054E+01 

0.9 

0.0 

-7.02472E+00 

-1 . 19310E+01 

2 

0.05000 

0.0 

-3.33523E+01 

-3.55047E+01 

-2.42968E+01 

0.0 

0.0 

-2.04410E+01 

-3.24692E+01 

1 

-0.05000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

X* 


0.0 

0.10000 

0.20000 

0.30000 

0.40000 

0.50000 

0.60000 

0.70000 

J 

V 

10 

11 

12 

13 

14 




14 

1.25000 

0.9 

0.0 

0.0 

0.0 

0.0 




13 

1.15000 

2 . 99216E+01 

1.73914E+0I 

0.0 

0.0 

0.0 




12 

1.05000 

2. 72814E+01 

1 .4464 IE +01 

0.0 

0.0 

0.0 




11 

0.95000 

3.452191+01 

1.60422E+01 

0.0 

0.0 

0.0 




10 

0.65000 

5.33483E+01 

2.26507E+01 

o.o 

0.0 

0.9 




9 

0.75000 

0.0 

6. 20654E+00 

0.0 

0.0 

0.0 




6 

0.45000 

2.64275E+00 

4.49058E+00 

0.0 

0.0 

9.0 




7 

0.55000 

5. 90830E+00 

6. 76106E+00 

0.0 

0.0 

0.0 




4 

0.45000 

1.32965E+01 

2.50802E+0! 

4.64201E+01 

2. 15346E+01 

2. 15499E+01 




5 

0.35000 

I.40931E+OI 

Z.54901E+01 

4. 25045E+01 

4. 06003E+01 

4. 06096E+01 




4 

0.25000 

4.63776E+00 

1 .68913E+01 

3. 13630E+01 

3. 73642E+01 

3. 73655E+01 




3 

0.15000 • 

■1 . 33304E+01 

2.46159E+00 

1.63610C+01 

2.64066E+01 

2.S3949E+01 




2 

0.05000 - 

•3.43213E+01 

-1 . 9646 IE +01 

1.21426E+01 

2.55433E+01 

2 . 55137E+01 




1 ■ 

•0.05000 

0.0 

0.0 

0.0 

0.0 

0.0 




X* 


0.60000 

0.90000 

1.00000 

1.10000 

1.20000 





Figure 20 

Complete Listing of U-Velocity Field Variable Calculated On All Planes 
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U VELOCITY 

J 

Y 

2 

3 

14 

1.25000 

0.0 

0.0 

13 

1.15000 

0.0 

1.159532*01 

12 

1.05000 

0.0 

1 ,08592E*01 

11 

0.95000 

0.0 

8.67430E*00 

10 

0.65000 

0.0 

5.620812*00 

9 

0.75000 

0.0 

2.38451E+00 

8 

0.65000 

0.0 

5. 36769E-01 

7 

0.55000 

0.0 

1.56484E+00 

6 

0.45000 

0.0 

7. 06035E+00 

5 

0.35000 

0.0 

1 . 36416E*0 1 

4 

0.25000 

0.0 

1.69720E+01 

3 

0.15000 

0.0 

1 . 66152E*01 

2 

0.05000 

0.0 

1.71299E*0i 

1 

-0.05000 

0.0 

0.0 

X* 


0.0 
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0.0 

0.0 

0.0 

2.155162*01 

4.061552*01 

3.73871E*01 

2.839772*01 

2.55156E*01 

0.0 

1.20000 


6 

0.0 

3.99235E*0l 

3.623112*01 

3.S0625E+01 

2.54897E+01 

0.0 

1.02796E+01 
8. 14635E+00 
6.851762*00 
4. 865122*00 
5.86673E-01 
-7. 01729E*00 
-2.06504E+01 
0.0 

0.60000 


9 

0.0 

3.702692*01 

3.38949E*01 

3.60360E*01 

3.440492*01 

0.0 

1.06595E*01 
1 . 029412*01 
1.234052*01 
1 . 05999E*0i 
1.871732*00 
-1 . 19336E*01 
-3.26633E*0I 
0.0 

0.70000 


J 

Y 


14 

1.25000 

0.0 

13 

1.15000 

0.0 

12 

1.05000 

0.0 

11 

0.95000 

0.0 

10 

0.85000 

0.0 

9 

0.75000 

0.0 

8 

0.65000 

0.0 

7 

0.55000 

0.0 

6 

0.45000 

0.0 

5 

0.35000 

0.0 

4 

0.25000 

0.0 

3 

0.15000 

0.0 

2 

0.05000 

0.0 

1 

-0.05000 

0.0 


AT Z( 8)* .65000 
6 


U VELOCITY 
3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.10000 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.20000 


5 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.30000 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.40000 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.50000 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.60000 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.70000 


Figure 20 (continued) 
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J Y 
14 1.25000 

13 1.15000 

12 1.05000 

11 0.95000 

10 0.05000 

9 0.75000 

0 0.05000 

7 0.55000 

0 0.45000 

3 0.35000 

4 0.25000 

3 0.15000 

2 0.05000 

1 > 0.05000 
X» 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.60000 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.90000 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

l. 00000 


0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1.10000 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1.20000 


Figure 20 (continued) 
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3.6 A Samp 1 e Terminal Sess i on 


TEACH can be run from any virtual machine that has 4.0 Megabytes of storage and is linked 
to disk CMB1— F by typing f, U57103 Tl . If a Tektronix terminal is used# the CalComp preview 
package di splays all plots on the screen making any CalComp plots optional. On any other 
terminal# all plots# except the convergence plot# are plotted by the CalComp. 

The menu that first appears will be a subset of the one shown below. Options will appear 
in the menu as they become relevant. For example# the M RC n option is not relevant until 
the input files have been identified (option "11") and is not displayed initially. 

A second filename and filetype (U57103 ENTHALPY in this case) is not required unless tem- 
peratures are being calculated. 




3-D TEACH 
SYSTEM MENU 


II 

IDENTIFY THE INPUT FILES 


El 

EDIT THE MAIN INPUT FILE 


RR 

READ RESTART FILE FROM PANVALET 


RC 

RUN ON CMS 


RB 

RUN ON BATCH 


E0 

EDIT THE OUTPUT 


PO 

PRINT THE OUTPUT 


FO 

FICHE THE OUTPUT 


R I 

RENAME THE MAIN INPUT FILE 


MI 

PUT INPUT CARDS ON MINIDISK 


PR 

PUT RESTART FILE ON PANVALET 


D5 

DISPLAY STATUS OF FILES BEING USED 


ST 

STOP 

CHOOSE AN 
RC 

OPTION: 


If the option to make temporary input revisions is chosen# a dummy file is created and the 
user is put into the edit mode to make modifications while the original file is left un- 
ci i sturbed . 

If the option to put restart file on PANVALET is chosen# a batch job is submitted. There- 
fore the file will not immediately be put on PANVALET. This delay should be taken into 
consideration when using this option so that files are not temporarily irretrievable. 


preceding page blank not F|LM£D 
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The menu appears on the screen whenever a new option is to be selected. Option 
selected first. If the same case had been run earlier and the results saved in 
file/ those results could be loaded back into the program by selecting option 4 
example option 3 was selected allowing the user to check his breakup before 
with the analysis. 


3-D TEACH PROGRAM 

> 1 - Read a case from cards 

3 - Display geometry breakup 

4 - Start From a restart file 

5 - Run case 

6 - Streak-lino plots 

8 - Sava fiald variables in restart file 

9 - Sava field variables in the output file 

13 - Call user supplied subroutine USER1 

14 - Call user supplied subroutine USERS 

99 - Stop 
Choose an option 


.3 

UHAT CELL PLANE (JOUID VOU LIKE TO UIEU? 

(IE. I«2# X*1 .5# K-3, Y-1.2E-1# STOP) 

.z - .3 

DO YOU UANT FULL SIZE PLOTS? .NO 

DO YOU UANT ONE-TO-ONE SCALING OF THE AXES? .NO 


11/16/84 

03/04/86 


1 must be 
a restart 
In this 
proceedi ng 
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It was determined that the geometry breakup is acceptable so option 5 is selected to run 
the case. 


3-D TEACH PROGRAfl 

* i - Read a case from cards 
> 3 - Display geometry breakup 

4 - Start from a restart file 

5 - Run case 

6 - Streak-line plots 

3 - Save field variables xn restart file 
9 - Save field variables in the output file 

13 - Call user supplied subroutine USERl 

14 - Call user supplied subroutine USERS 

99 - Stop 
Choose an option 


• 5 
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The convergence plot is drawn as the case proceeds showing the values of each residual 
source as they approach the goal (SORMAX). If after several iterations it does not look 
as though the case will converge the iteration loop may be exited by hitting^ carriage re 
turn to put the terminal into CMS and typing "HALT”. In this case computation will stop 
after the next set of residual sources are plotted. In any case, after the iterations are 
complete, the menu will reappear. 



SfSSSS 

SSSffS 

, '<14 
. T4P 
♦ OIV 
. 'I Lie 

. tmz 

, OCmS 
. *#«* 

. TIM 


1.0 TTCM 
oa .14 


TYPE 'C' FOR A CALCOflP COPY - ELSE HIT RETURN TO CONTINUE. 
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3-D TEACH PROGRAM. 

t 1 - Read a case from cards 
* 3 - Display geometry breakup 
4 - Start From a restart file 
> 5 - Run case 

6 - Streak-line plots 

8 - Save field variables in restart file 

9 - Save field variables in the output file 

13 - Call 
H - Call 

99 - Stop 

Choose an option 

user supplied subroutine USER1 
user supplied subroutine USERS 

.6 



The streakline plot is an approximate graphic representat i on of what might be seen in a 
water tunnel when particles in the fluid are illuminated as they pass through a fixed 
plane. A plane perpendi cular to any of the coordinate axes and located anywhere in the 
flow field is selected by the user. One hundred massless particles are injected into the 
flow field from random locations in the plane and their path is approximated using the 
Runge— Kutta predi ct i on— cor rector method. The number of time steps and the time interval 
of each step can be selected by the user. If not, the number of time steps defaults to 
twenty and the time i nterval defaults to a value based on the highest velocity in the 
field and diagonal length of the control volume. Care should be exercised in choosing the 
number of time steps and the time interval. Too many time steps and/or too large a time 
interval will have the effect of making the ’laser beam’ quite thick and the streak line 
plot will not be representat i vg of the flow in a plane but rather flow in a block. One 
has to be especially careful at the exit. In interpreting the st reak 1 i nes i t should be 
kept in mind, that at best, this is a qualitative picture. Quantitative information 
should never be extracted from this plot. 
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Option 9 places all field variables in the output file for printing later. 

Option S is used to save all the field variables calculated by the program for restarting 
at a later tine. This allows the user to continue running a case that had not _ converged 
yet, to plot the results in a later session, or to change the boundary conditions and 
start using the restart file as initial guesses thereby saving CPU time. Option 99 takes 
the user out of the 3D— TEACH program menu and returns to the system menu. In this menu 
the saved output file can be edited and printed if need be by using options E0 and P0 re- 
spectively. Sometimes the output file becomes too large to be edited inside the menu. In 
that case, it can usually be edited by coming out of the menu by typing ST and typing XE 
U571 OUTPUT D. 



3-D TEACH PROGRAM 


t 1 - Read a case from cards 
I 3 - Display geometry breakup 
A - Start from a restart fi\e 
* 5 - Run case 

> S - Streak-Une plot* . . . 

g _ g ave field variables in restart file 
9 - Save field variables in the output file 

13 - Call user supplied subroutine USER1 

14 - Call user supplied subroutine USERS 

99 - Stop 
Choose an option 


. 8 

Results saved in restart file 
Choose an option 


.9 

Output has been saved on disk 
Choose an option 


99 
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Occasionally a user is interested in printing or plotting parameters which are not calcu- 
lated by the production version of TEACH. Two dummy subroutines in the prodcution 
version, called USER1 and USER2 , have been provided as "place holders" to give 
cess to all the TEACH variables for calculation of special parameters. Either or 
these subroutines can be replaced by subroutines written by the user and 
lecting option 13 or 14 from the MAIN Menu ("Call user supplied subroutine 
user supplied subroutine USER2" respectively). Because individual needs 
planat i on of how to write a user supplied subroutine will be given in the 
of a programmer familiar with TEACH should be solicited. 


users ac- 
both of 
i nvoked by se- 
USER1" or "call 
may vary, no ex- 
manual . The ai d 


3-D TEACH PROGRAM.. 

1 - Read a case from cards 

3 - Display geometry breakun 

4 - Start from a restart file 

5 - Run case 

6 - Streak-line plots 

8 - Save field variables in restart file 

9 - Save field variables in the output file 

13 - Call user supplied subroutine USER1 

14 - Call user supplied subroutine USERS 

99 - Stop 
Choose an option 


.11 
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4 . 0 SAMPLE CASES 


Three cases are given in this section to give the user an idea of the range of problems 
that can be solved by TEACH. Much effort has been expended to make the input structured 
and logical. After setting up a few cases the user should be able to set up and modify 
new cases with little help from the manual. The three cases are arranged in order of com- 
plexity starting with a very simple case (Laminar Flow in a Pipe, case 17) to one having 
complex boundary conditions (JT9D Louver, case 23). The user should go through each case 
in order before attempting to set up a case on his own. Each sample case is chosen be- 
cause of its similarities to and differences from the previous case so that an 
understanding of the philosophy behind the input can be gained a step at a time. 

4 . 1 Laminar Flow in a Pipe 


4.1.1 P rob 1 em 


To model laminar airflow in a .1 ft diameter by 3.2 ft long pipe. The air at 15 psia and 
70°F is entering the pipe with a uniform velocity profile of .1 ft/sec. 

4.1.2 Input 


This is a simple problem which can be modeled with one block. As a first pass, 20 cells 
are desired in the X— direction and 8 cells in the Y— direction. Because the velocity will 
be changing more rapidly in the vicinity of the inlet and more slowly downstream of the 
inlet, a finer grid has to be used in the inlet region. An expansion factor of 1.1 will 
provide this adequately (card 3A). As the boundary layer builds up with increasing dis- 
tance from the inlet, some high velocity gradients will be present near the wall, 
requiring a fine breakup near the pipe wall. An expansion factor of .9 in the Y— direction 
will give the extra resolution needed near the wall (card 3B). There are no variations in 
properties in the Z— direction so making the section of the pipe being modeled three radi- 
ans thick is arbitrary. One block has five boundary condition card pairs so five boundary 
condition card pairs (cards 5A and 5B) are needed. The first card pair describes the ve^ 
locity, temperature, pressure, and the coefficient of turbulent kinetic energy for the 
inlet. The second card pair makes the east wall of the block an outlet in which fluid 
properties must be calculated by the program. The next card pair places a wall having a 
temperature of 70°F at the north boundary, the fourth pair makes the south boundary an ax- 
is of symmetry, and the last pair makes thQ front and back walls periodic. 

The case indicator card (card 6) is set up to treat the problem as a body of revolution. 
All the field variables will be initialized to the value of cell (1,4,2) in the inlet and 
cell (20,6,4) will be monitored in the output. None of the other variables on the card 
require values for the case. Case indicator card 7 states that the flow is laminar, den- 
sity is constant, the viscosity to be used is input, all input and output units are Eng- 
lish, and temperature is constant. The card also sets a maximum residual source of .005 
and maximum iterations to 100. 

The fluid properties card (card 8) gives some of the constants required by the program. 
The calculation switches card (card 9) has turned on U—veiocity, V— velocity, W— velocity, 
and pressure calculations. All the other cards in the input are left blank allowing tho 
default values to be set by the program. 


ORIGINAL PAGE IS 
OF POOR QUALITY 



TEACH 3D 


PAGE 83 


ISSUED: 11/16/84 
REVISED: 03 / 04/86 



Figure 21 

Input for Laminar Flow in a Pipe 
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4.1.3 Output 


Figures 21 and 22 show the card input before 
respectively. Figure 23 shows the geometry 
TEACH can approach, with a sufficiently fine 
tents and purposes, be considered exact. In 
results. The exact solution is a parabolic v 
walls and a velocity of twice the average on 


and after the program has added the defaults, 
plot. This case has an exact solution which 
breakup, to a point which could, for all i n- 
thi s case, even a coarse grid gives very good 
elocity field having a velocity of 0.0 at the 
the centerline (Figure 24). 


INPUT WITH DEFAULTS: 


EXAMPLE 4. 

1 LAMINAR 

FLOW IN A 

PIPE 





1 

i 1 

l 1 

5 10 

0.0 

0.0 

0.0 

0 0 


2 

3.2000 

20 1.10 







3A 

0.0500 

8 0.90 







3B 

3.0000 

3 1.00 







3C 

1 1 

1 1 

2 4 

3 5 

5 




4 

1 SOP 

0.10 

0.0 

0.0 

15.00 

70.00 

0.01 

0 1.00 

5A 

0.0 

0.0 

0.0 

0.0 




48000. 

5B 

2 UOP 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 1.00 

5A 

0.0 

0.0 

0.0 

0.0 




0 . 

5B 

3 SWL 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 1.00 

5A 

0.0 

0.0 

0.0 

0.0 




0 . 

5B 

4 SYM 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 0.0 

5A 

0.0 

0.0 

0.0 

0.0 




0. 

5B 

5 PER 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 1.00 

5A 

0.0 

0.0 

0.0 

0.0 




0. 

5B 

2 1 

4 2 

20 6 

2 

100. 

100. 

1 . 

1000000. 

6 

T 0 

T 1000 

2 500 

0.0050 

-30 0 

0 

-1.00000 

0.0 

7 

0.0050 

. 1110E-G4 , 

.8050E-01 

29.00 

0.700 

1.44000 

0.0 

0.25000 

8 

T T 

T T 

F F 

F F 

F F 

F 



9 

0.5000.5000 

.5001.0000. 

.7000.7001. 

0000.7000 

.3000.7000 

.8001.000 



10 

3 3 

3 5 

3 3 

3 3 

10 




11 

0.0900 

1.0000 

1.4400 

1.9200 

0.4187 

9.7930 

1.0000 

0.9000 

12 









13A 

F 

F 

F 

F 




1 

13B 

2.01600 

32.00000 

44.00999 

18.01601 





13C 

I. 00000 

l. 00000 

1.00000 

1.00000 





13D 

0.90000 

0.90000 

0.90000 

0.90000 





13E 

1 2 

3 4 


7.9370 

1.0000 

4.0000 

2.0000 


14 


FIGURE 22 

Input with Defaults for Laminar Flow in a Pipe 
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Z=l -S 

K( CELL )= 3 
K( BLOCK ) = l 


'-b.22 o’.ia 0.64 0.92 


u5t 7Ti z ’.os FITS 2 ’. so 5775 F.w 

X (FT) 


FIGURE 23 

Grid Breakup for Laminar Flow in a Pipe 
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ommm, page ic 

OF POOR QUALITY 


J Y 
10 0.05109 

9 0. 04601 

8 0,04359 

7 0.03366 

6 0 . 0 3322 

5 0.02716 

4 0.02042 

3 0.01294 

2 0.00462 

1 - 0.00462 
Xs 

J Y 
10 0,05199 

9 0,04801 

8 0.04359 

7 0.03666 

6 0.03322 

5 0.02716 

4 0.02042 

3 0.01294 

2 0.00462 

1 - 0.00462 
X* 

J Y 
10 0.05199 

9 0.04601 

8 0.04359 

7 0.03866 

6 0.03322 

5 0.02716 

4 0.02042 

3 0.01294 

2 0.00462 
1 - 0.00462 

X* 


2 

0.0 

1 . 00000E-01 
1 .OOOOOE-Ol 
1.000006-01 
1 . 000006-01 
1 . OOOOOE-Ol 
1. OOOOOE-Ol 
1. OOOOOE-Ol 
1 .OOOOOE-Ol 
1. OOOOOE-Ol 
0.0 
10 

0.0 

1 . 59520E-02 
4 . 83948E-02 
6 .072 1 OE-02 
I . I2027E-01 
1.41 I02E-01 
1 .66371E-01 
1 . 85859E-01 
1 . 971 37E-01 
1.97137E-01 
0.63894 
18 

0.0 

1.58931E-02 
4.82439E-02 
8.05262E-02 
1. 11874E-01 
1.411 24E-01 
1 . 66734E-0 1 
1.86684E-01 
1 . 98351E-01 
1 . 96352E-0I 
2.00655 


U VELOCITY 
3 

0.0 

2 . 78490E-02 
6.97924E-02 
9.95799E-02 
1.19471E-01 
1 . 31826E-01 
1.35823E-01 
1.42305E-01 
I.43669E-01 
1.43669E-0I 
0.05567 
11 

0.0 

1 . 59200E-02 
4.83I09E-02 
6.0609IE-02 
1.11931E-01 
1.41069E-01 
1 . 66525E-0 1 
1 . 86243E-01 
1 . 97714E-0I 
1.97714E-01 
0.75870 
19 

0.0 

1 . 56895E-02 
4.62339E-02 
6.05107E-02 
1 . 1 1654E-01 
1 .41 100E-01 
1 . 66708E-01 
l . 86657E -01 
1.98324E-01 
1.98324E-01 
2.26526 


4 

0.0 

1 . 929026-02 
5.66768E-02 
9.079786-02 
I. 16696E-01 
l . 39513E-0I 
1 . 53521E-0I 
1.61621E-01 
l .65664E-01 
1 .65664E-01 
0.11733 
12 

0.0 

1 . 590466-02 
4.62716E-0Z 
6. 05575E-02 
1 . H669E-01 
1.41090E-01 
1 .666141-01 
1 .66451 E-0 l 
1.98023E-01 
1.960Z3E-01 
0,69044 
20 

0.0 

1.5664 IE-02 
4.62 I62C-02 
8.04658E-02 
1 . 1 1621 E-0 1 
1.41061E-01 
1 . 66663E-0I 
1.86609E-01 
I .96274E-01 
1.96275E-01 
2.54767 


5 

0.0 

1 . 74976C-0Z 
5.24394E-02 
6. 59612E-02 
1.16157E-01 
1.411736-01 
1.59662E-01 
1.7Z109E-01 
1. 78233E-01 
1 . 78234E-01 
0,18493 
13 

0.0 

1 .569851-02 
4.62556E-02 
8.053796-02 
1 . U676E-01 
1.41 100E-01 
1 .666666-01 
i .64564E-0I 
1 . 96166E-01 
1.96167E-01 
1.03536 
21 

0.0 

1.56766E-02 
4.61966E-02 
6. 04508E-02 
1 . J 1 774E-0I 
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FIGURE 24 


Axial Velocity Field for Laminar Pipe Flow 
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A . 2 Turbulent FIom in a Channel 

4 . 2 . 1 Problem 


To model turbulent flow in a .05 ft wide by 4. ft long channel. The air at 15 psia and 
7 0 0 F is entering the channel with a uniform velocity profile of 93. ft/sec and an energy 
dissipation of 48261 ft 2 /sec 2 . The air at the inlet has a Re = 67446, which means that 
the flow is turbulent. 


4.2.2 Input 


The block breakup is the same as the previous sample case with the exception that a finer 
grid (25X16) is being used and there is no axis of symmetry, so the boundary condition de- 
scription cards for that are eliminated. The north and south wall use the same boundary 
condition description card pair (cards 5A and 5B). The case indicator cards have been 
changed to indicate turbulent flow, with viscosity calculated every iteration. The maxi- 
mum number of iterations allowed has been increased to 200 because of the finer breakup 
being used. The fluid properties (card 8) are the same as the previous case, so no change 
is made to this card. The calculation switches card (card 9) has been changed to turn on 
the calculation of kinetic energy and energy dissipation. The rest of the input is left 
blank allowing default values to be set by the program. 

Figures 25 and 26 show the input without and then with defaults. 


4.2.3 Output 


Figure 27 shows the geometry breakup. All goemetry breakup plots show the boundaries of 
the scalar cells which are bounded on the west and east by u-velocity gridlines and on the 
south and north by v-velocity gridlines. Figures 28 and 29 show the U-velocity and pres- 
sure fields respectively (only the values near the inlet and exit are reproduced here). 
Note that the maximum velocity (103 ft/sec) is only 102. higher than the the average veloc- 
ity (93 ft/sec) when the flow is turbulent. 
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FIGURE 26 

Input with Defaults for Turbulent Flow in a Channel 
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9.166182*01 9.169432*01 9.190902*01 
6.654042*01 8.657242*01 6.65816E*01 
7.921862*01 7.924372*01 7.924532*01 
6.645992*01 6.647562*01 6.64710E*0I 
0.0 0.0 0.0 

2.60314 2.90413 3.23522 


Axial Velocity Field for Turbulent Channel Flow 


TEACH 3D 


PAGE 92 


ISSUED: 11/16/8A 
REVISED: 03 / 04/86 

ORIGINAL PAGE IS 

OF POOR QUALITY 


STATIC PGCSSLWC 


( LB/IN* ■ 2 ) 


AT 21 3> = 1.5000 


J Y 
! J o.c5io4 
17 o.o- 3?* 
0 . 0466 7 
15 ?.C4415 

iA o.:;n 8 

! 3 0.03533 

3.034*3 

u o.::i*9 
: 3 a . ::;:3 
’ c.o::77 

6 0.61071 

7 0.01502 

i C.CII67 
5 0.C0C62 

4 0.005*5 

3 0.00333 

: C.00U4 

1 -C. 30104 


1 

0.0 

1 .500012*01 
1 . 50C0 1C *01 
1.500012*01 

i.soeoic«oi 

1 . 50001 C *01 
I .500015*01 
1.500012*01 
l .500015*01 
1 .5000 12*01 
1 .5000 12*01 
1 .500012*01 
1 . 5C00 1 C *0 1 
1 - 50001 C«OK 
1 .500012*01 
1 . 500012*01 
1.560012*01 
0.0 

-0.01937 


2 

0.0 

1.496172*01 
1 .496232*01 
1.496232*0! 
1 .496222*01 
1.496222*01 
1.493222*01 
1.493212*01 
1.496212*01 
1.496212*01 
1.499212*01 
1.496222*01 
1.496222*01 
1.493222*01 
1.496232*01 
1.493232*01 
1.496172*01 
0.0 

0.01937 


3 

0.0 

1.496062*01 
1 .496122*01 
1.496122*01 
1.498122*01 
1.496122*01 
1.496122*01 
1.496122*01 
1.496122*01 
1.496122*01 
1.496122*01 
1.496122*01 
1.496122*01 
1.496122*01 
1.496122*01 
1.496122*01 
1 .496062*01 
0.0 

0.04196 


4 

0.0 

1.497942*01 
1.496012*01 
1.493012*01 
1.496012*01 
I .498012*01 
1.498012*01 
1.498012*01 
1.498012*01 
1.496012*0! 
1.496012*01 
1.496012*01 
1.496012*01 
1.496012*01 
1.493012*01 
1.493012*01 
1.497942*01 
0.0 

0.10665 


5 

0.0 

1.497632*0! 
1.497652*01 
1 .497692*0! 
1.497692*01 
1 .497692*01 
l .497632*0 I 
1.497692*01 
1 .497692*01 
1.497692*01 
I .497392*01 
1 .497392*01 
1.497392*01 
1.497392*01 
1. 497692*01 
1.497662*01 
1.497632*01 
0.0 

0. 16040 


6 

0.0 

1.497702*01 
1.497752*01 
1.497772*01 
1,497772*01 
1,497772*01 
1,497772*01 
l .497762*01 
1.497762*01 
1.497762*01 
1.497762*01 
1.497772*01 
1.497772*01 
1 .497772*01 
1.497772*01 
1.497752*01 
1.497702*01 
0.0 

0.71712 


7 

0.0 

1.497572*01 
1.497612*01 
1.497632*01 
1.497632*01 
1 .497635*01 
1.497632*01 
I ,497632*01 
1.4976:2*01 
1.497632*01 
1.497632*01 
1.497632*01 
1.497632*01 
1.497632*01 
1 .497632*01 
1.497612*01 
1.497572*01 
0.9 

0.27950 


6 

0.0 

1 .497432*01 
1 .497462*01 
1 .497492*01 
1 .497452*01 
1.497502*01 
1.497502*01 
1.497502*01 
1.497502*01 
1.497502*01 
1.497502*01 
1.497502*01 
1 .497502*01 
I .497492*01 
1.497492*01 
1.497462*01 
1.497432*01 
0.0 

0.34612 


J Y 

15 0.05104 

17 0.04396 

16 0.04667 

15 0.04415 

14 0.04136 

13 0.03333 

12 0.03498 

11 0.03129 

10 0.02723 

9 0.02277 

3 0.01871 

7 0.01502 

6 0.01167 

5 0.00862 

4 0.00535 

3 0.00333 

2 0.00104 

1 -0.00104 


25 

0.0 

1.49255E+01 
1.49256E+01 
1.49257E+01 
1.49258E401 
1.49259E+01 
1 .49260E+01 
1 . 49261E+01 
1 .49262E+0I 
1 . 49262E+01 
1.49261E+01 
1 .49260E+01 
1.49259E+0I 
1.4925QE+01 
1 . 49257E+01 
1 . 49256E+01 
1 .49255E+01 
0.0 

3.40864 


26 

0.0 

1.49195E*01 
1.49196E+01 
I.49197E+0I 
1.49198E+01 
1.49199E+01 
1 .49200E+01 
I.49201E+01 
1 .49202E+0 1 
I.49202E+01 
1 . 4920 1 E+O 1 
1.49200E+01 
1 .49199E+01 
1 . 49I98E+01 
1 . 49197E+01 
1 .491 96 E+ 01 
1 .49195E+01 
0.0 

3.79018 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

4.20987 


FIGURE 29 


Pressure Field for Turbulent Channel Flow 
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4.3.2 Input 

Figure 31 gives the cell breakup in the plane perpendicular to the flow. In order to ac- 
curately calculate the boundary layers, a finer breakup was used near the south and north 
walls (y=-. 03125 ft and y - . 03125 ft). Uniform breakups were used in the X and Z di- 
rections. Cross-sections of other directions are not shown because of the high aspect ra- 
tio of the length to height and length to width. Figure 32 shows the input for this case* 


I 



“0.02 -p.QO O.Ol 0.03 0-05 ^ 0.06^ 0*08 0.09 
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J 

2345 6 7 8 9 1 011 12 1 3 1 4 151611 - 8 - 9 -^ 

mil 1 1 1 1 1 1 2 2 2 2 2 2-2 ^> <> c 

<< + 

1 w 

If II 1' I I I I I I I — I 1 1 — i — i — i — m 


1=2 

I ( CELL )= 2 
I ( BLOCK )= 1 


ISSUED: 

REVISED: 


0.04 


11/16/84 

03/04/86 
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CARD INPUT: 


3-0 TEACH 

EXPERIMENTAL DUCTU. 

MOORE) / AR 

= 4. / BELOW PUNE 

OF SYMMETRY 

1 

2 z 

1 4 

6 10 

-3.00000 

-0.03125 

0.0 

0 1 


2 

0.0 

6 1.00 







3A 

3.0000 

6 1.00 







3A 

0.0 

il 1.15 







3B 

0.0312 

7 0,83 







3B 

0.1250 

6 1.00 







3C 

1 1 

1 16 

0 1 

0 1 

5 




4 

2 1 

l 0 

13 1 

0 1 

5 




4 

1 2 

1 16 

0 0 

1 1 

5 




4 

2 2 

l 0 

13 0 

1 1 

5 




4 

1 SUL 

0.0 

0.0 

0.0 

0.0 

167.00 

0.01 

0 0.0 

5A 

0.0 

0.0 

0.0 

0.0 




0 . 

5B 

5 SYH 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 0.0 

5A 

0.0 

0.0 

0.0 

0.0 




0 . 

5B 

13 UOP 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 0.0 

5A 

0.0 

0.0 

0.0 

0.0 




0 . 

5B 

16 SOP 

55.00 

0.0 

0.0 

14.70 

60.00 

0.01 

0 0.0 

5A 

0.0 

0.0 

0.0 

0.0 




453. 

5B 

1 1 

2 2 

0 0 

0 

100. 

100. 

l. 

1000000. 

6 

F 0 

F 0 

2 800 

0.0200 

-10 0 

0 

0.0 

0.10 

7 

i .0 

.0 

18.21 





1 

1 7C 

0.0050 

. 1200E-04 

.7700E-01 

29.00 

0.700 

i. 00000 

0.0 

0.24000 

8 

0T -IT 

-IT -IT 

-IT -IT 

OF OF 

OF OF 

OF 



9 

0.0 0.0 0 

.0 0.0 0 

.0 0.0 0 

.0 0.0 0. 

0 0.0 0. 

0 0.0 



10 

0 0 

0 0 

0 0 

0 0 

0 




11 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

12 









13A 

F 

F 

F 

F 




0 

13B 

0.0 

0.0 

0.0 

0.0 





13C 

0.0 

0.0 

0.0 

0.0 





13D 

0.0 

0.0 

0.0 

0.0 





13E 


FIGURE 32 


Input File for John Moore’s Experiment 
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3 . 3 Output 


gure 33 shows the axial velocities in the duct at the plane 
column represent the velocities as they vary from the south 
use of Co r i o 1 1 i s forces acting on the gas. The velocity is 
en by locating the highest velocity at each axial station. 


of symmetry 
wall to the 
peaking off 


The numbers i n 
north wa 11. B 
center as can 


J r 
20 0.03218 

19 0.03012 

18 0.02740 

17 0.02412 

16 0.02017 

15 0.01542 

14 0.00966 

13 0.00278 

12 -0.00278 

11 -0.00761 
10 -0.01161 

9 -0.01546 

8 -0.01663 
7 -0.02140 

6 -0.02380 

5 -0.02589 

4 -0.02770 
3 -0.02928 

2 -0.03065 

1 -0.03185 
Xs 

J T 

20 0.03238 

19 0.03012 

16 0.02740 

17 0.02412 

16 0.02017 

15 0.01542 

14 0.00966 

13 0.00278 

12 -0.00278 
11 -0.00761 
10 -0.01181 

9 -0.01546 

3 -0.01863 

7 -0.02140 

6 -0.02380 

5 -0.02589 

4 -0.02770 
3 -0.02928 

2 -0.03065 
1 -0.03165 

X= 


2 

0.0 

5.500002*01 
5 . 500002*0 1 
5.500002*01 
5 . 50000E +01 
5 . 500002*01 
5 . 50000E *01 
5 . 50000E *01 
5.500002*01 
5 . 50000E+0 1 
5 . 500002*01 
5.500002*01 
5 . 50000E+0! 
5 . 500002*01 
5.500002*01 
5 . 50000E*0 1 
5.500002*01 
5.500002*01 
5.500002*01 
0.0 

- 3.00000 

10 

0.0 

4.107182*01 
4 . 86502E+01 
5.316732*01 
5.6301 1E*0! 
5.881972*01 
6 . 08812E*0 1 
6.239822*01 
6 . 29334E*0 l 
6 . 28531 E *0 1 
6 . 22670E *0 1 
6. 12534C*01 
5.988972*01 
5.822432+01 
5.625202*01 
5.383662+01 
5.09035E+01 
4.670452*01 
3. 8996 1 E*0 l 
0.0 

l .00000 


U VELOCITY 

3 

0.0 

4.3201 9E+01 
5.170812*01 
5.52985E+01 
5.681282+01 
5.753242*01 
5.78945E+01 
5. 804472*01 
5. 80341E+01 
5.793172*01 
5.77452E+01 
5.746102*01 
5.704572*01 
5.643732*01 
5 . 551 99E*0l 
5. 40665E *01 
5. 159462*01 
4 . 70677E *01 
3.83904E*01 
0.0 

-2.50000 

11 

0.0 

4.105152*01 
4.881932*01 
5.312872*01 
5.625472*01 
5.877092*01 
6 . 08443E+01 
6.240712*01 
6.300632*01 
6.299922*01 
6.248382*01 
6.152482*01 
6.019302*01 
5.854112*01 
5.656832*01 
5.419402*01 
5. 119492*01 
4.697142*01 
3.921 75E*0 1 
0,0 

1.50000 


0.0 

4.131952*01 
4.972932*01 
5.427882*01 
5. 70555E+01 
5. 86898E+01 
5.958362*01 
5 . 99751E*01 
5.99758E*0l 
5.9757lE*0l 
5.933512*01 
5.868072*01 
5.773912*01 
5.643632*01 
5.469372*01 
5.239762*01 
4.936382*01 
4.51052E*01 
3.751392*01 
0.0 

- 2.00000 

12 

0.0 

4.105622*01 
4.662172*01 
5.312732*01 
5.624962*01 
5.67648E*01 
6.064662*01 
6.244072*01 
6.308472*01 
6.313102*01 
6.266942*01 
6. 175312*01 
6.044662*01 
5 . A8047E+01 
5.663102*01 
5.444672*01 
5.143592*01 
4.719222*01 
3.940132*01 
0.0 

2.00000 


0,0 

4.116532*01 
4.916732*01 
5.366732*01 
5.682852*01 
5 . 9071 1 E*01 
6.051062*01 
6.12156E+0I 
6.12S63E*01 
6.092272*01 
6.026292*01 
5.928032*01 
5.797662*01 
5.635452*01 
5.439672*01 
5.203992*01 
4.908222*01 
4.497052*01 
3.751252*01 
0.0 

- 1.50000 

13 

0.0 

4.106562*01 
4.665462*01 
5. 31626E*0L 
5.626502*01 
5.880132*01 
6.088942*01 
6.250362*01 
6.317632*01 
6.325712*01 
6.28310E+01 
6.194422*01 
6.065572*01 
5.902042*01 
5.704552*01 
5.465672*01 
5.163272*01 
4.737112*01 
3.954722*01 
0.0 

2.50006 


AT Zl 81* 

6 

0.0 

4.115682*01 
4. 90435E+01 
5.346352*01 
5.666342*01 
5.912652*01 
6.041272*01 
6.191542*01 
6.203672*01 
6.165132*01 
6.065892*01 
5.972022*01 
5.626812*01 
5.659122*01 
5.461332*01 
$.226962*01 
4.934192*01 
4.524692*01 
3.776872*01 
6.0 

- 1.00000 

14 

0.0 

4.118802*01 

4.695612*01 

5.326592*01 

5.638842*01 

5.890472*01 

6.699292*01 

6.260742*01 

6.328002*01 

6.336072*01 

6.29344E*01 

6.204752*01 

6,075902*01 

5.91237E+01 

5.714902*01 

5.476022*01 

$.173622*01 

4.747462*01 

3.96$062*01 

6.6 

3.00600 


.13542 

7 

0.0 

4.116202*01 

4.900562*01 

5.338102*01 

5.656212*01 

5.908022*01 

6.102632*01 

6.224442*01 

6.247602*01 

6.212232*01 

6.131572*01 

6.014992*01 

$.869872*01 

5.699512*01 

$.501902*01 

5.267832*01 

4.974712*01 

4.563592*01 

3.810142*01 

0.0 

-0.50000 


8 

0.0 

4.114522*01 

4.89620E*0! 

5.331062*01 

5.647102*01 

5.899622*01 

6.100892*01 

6.236722*01 

6.271422*01 

6.244442*01 

6.169352*01 

6.056182*01 

5.913022*01 

5.743662*01 

5.544242*01 

5.311582*01 

5.016922*01 

4.602912*01 

3.843312*01 

0.0 

0.0 


9 

0.0 

4.110892*01 

4.890422*01 

5.323462*01 

5.63791E+01 

5.89014E*01 

6.094502*01 

6.23967E*01 

6.284742*01 

6.267622*01 

6.200742*01 

6.093292+01 

5.953402+01 

5.785442+01 

5.588122+01 

5.35255E+01 

5.056042+01 

4.638982+01 

3.873442*01 

0.0 

0.50000 


FIGURE 33 


Axial Velocities in the Plane of Symmetry 


cr re 
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Var i ables 

Coefficient of the finite difference equation 
Area of a cell face 
Convection coefficient 
Diffusion coefficient 

Weighting factor in the finite difference equation 
Weighting factor in the finite difference equation 
Total energy 

i Specific enthalpy of the species 

Thermal conductivity 
1 Mass source 

i Mass fraction 

’ Static pressure 

\ Universal gas constant 

i Source term in the differential equation 

Temperature 
: Time 

i Velocity in the x-direction 

t Velocity vector 

f Velocity in the y-direction 

4 Velocity in the w-direction 

< Ax i al di stance 

f Radi al di stance 

z Tangential distance 

T Exchange coefficient 

x Kinematic viscosity 

3 Densi ty 

T Stress 


11/16/84 

03/04/86 


0 - 2 - 
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B 

b 

E 

6 

F 

f 

eff 

i 

j 

k 

1 

N 

n 

P 

S 

5 

w 

w 

X 

y 

z 


Subscripts 


Back node 
Back face 
East node 
East face 
Front node 
Front face 
Effect i ve 

Cell X-direction index 

Cell Y-direction index 

Cell Z-direction index 

Spec i e 

North node 

North face 

Node being computed 

South node 

South face 

West node 

West face 

Ax i al di rect i on 

Radial direction 

Tangential direction 


number 

number 

number 


11/16/84 

03/04/86 
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igndix B - Graph Paner Codes 


?de 

No . 

Descr i pt i on 

Width 

L 

605 

10 x 10 to the centimeter 

11" 

? 

602 

10 x 10 to the half inch 

11" 

5 

500 

Plain White 

34" 

t 

505 

10 x 10 to the centimeter 

34" 

5 

502 

10 x 10 to the half inch 

34" 

s 

501 

10 x 10 to the inch 

34" 

7 

GC 1381 3 

Probability Paper-perf orated 

11" 

5 

6 05B 

Centimeter - perforated 

11" 

9 

500 

My lar 

34" 

0 

600 

Plain White 

11" 

1 

GC 1 38 0 9- 1 08 1 7 

A size half inch - perforated 

11" 

*9 

Cm 

600 

Mylar 

11" 

3 

GC1 3812-1 0818 

A size centimeter - perforated 

11" 

4 

622 

2-cycle full log 

11" 

5 

640 

4-log cycles x 10 div. 

11" 

6 

GC 1 38 07 

3 cycles - perforated 

11" 

7 

GC1350S 

5 cycles - perforated 

11" 

3 

601 

10 x 10 to the inch 

11" 

9 

600 

With ballpoint pens 

11" 

1 

10819 

Weibull Prob. - perforated 

11" 

2 

600 

Mylar with Ballpoint pens 

11" 

9 

500 

Mylar with Ballpoint pens 

34" 

0 

500 

With ballpoint pens 

34" 

3 

500 

Vel lum34" 



OTE: Should this parameter be omitted or invalid, graph paper code Number 02 (10 x 10 to 
the half inch, 11" wide) will be used when making the plot. 
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6 . 3 Appendix C ~ Flow Diagram For Calculating Potentially Unstable Problems 


For complex problems the computer program may be unstable if all the flow complexities are 
T turned on 1 at once. The following guidelines are meant to guide the user in obtaining a 
stable solution. The Flow Diagram provided below is only a guide and does not guarantee a 
stable solution. 
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4 Appendix D Warnings and Rgmgdi P5 


he following list of warnings apd/or corresponding remedies have been collected as a re- 
ult of extensive use of 3D TEACH. Their appearance here is to help the user with common 
ccuring user solvable problems. 

) Because initial guesses for enthalpies may not be consistent with the initial guesses 
or specie mass fractions# temperatures and densities calculated based on these could 
ause divergence of the solution procedure. To avoid this problem# density and temper 
ture calculations are not turned on for at least 10 iterations. 

) If# during the iterative process# a divide check occurs in LI50LV# it is possible that 
11 flow i 5 going into or coming out of a single cell. If this is the case the user 
hould try perturbing the solution process by one of the following: a) put a contraction 
t the exit b) change the number of sweeps or the under-relaxat i on factor on the velocity 
calculations c) temporarily change the velocities on specified boundaries. 

,) Any blockage at an inlet must be at least two cells thick. 

) An unspecified opening can not be used in conjunction with an adjoining blockage. 

i) Ho flow should be entering the calculation domain from an unspecified opening and the 
: 1 o w should be parabolic. If negative flow is generated block part of the opening to ac- 
:elerate the flow and make all flow positive. 

Stability and accuracy of the flow field depends on cell aspect ratios and the expan" 
5 i on ratio. Cell a spect- rat i o s greater than 10.0 will give rise to instabilities. Any 
/alue greater than 1.0 will give rise to inaccuracy. This error will increase as the as- 
pect ratio is increased. Expansion ratios greater than 1.25 and less than 0.8 should also 
:e a vo i ded . 

7 ) When using streakline plots it should be kept in mind that the combination of the size 
>f the time step and the number of time steps should not be very large. A rule of thumb 
to follow would be that this time should not be more than one tenth of the time required 
for a particle to travel from one end of the geometry to the other. Special care should 
3 e taken when plotting streaklines near an outlet. In that case the time step should be 
Less than the time required for the particle to exit the geometry. 
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6 . 5 Appendix E ~ Program Maximums 


V a r i a b 1 e 

Cord i 

i t appears on 

Current Maximum Value 

NXBIK 


( Card 

Type 

2) 

58 

NYBLK 


(Card 

Type 

2) 

58 

NZBLK 


( Card 

Type 

2) 

16 

NFBLK 


(Card 

Type 

2) 

1500 

NBC 


(Card 

Type 

2) 

120 

Sum of 

NX 

( Card 

Type 

3A ) 

58 

Sum of 

NY 

(Card 

Type 

3B ) 

58 

Sum of 

NZ 

(Card 

Type 

30 

16 
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